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tDCK-nttKn.u 


ji semi-empiricaC unified mobiCity modeC dedicated to the modeCing of short channeC Cength 
MOS ( F c E < Ts for Cow power anaCog appCications has 6een deveCoped in this wor^. tfhe existing <BSI!M 
mobiCity modeC ps tafien up and modified here in order to incorporate the Coutomh scattering effect , 
which pCays the dominant roCe in determining the effective mobiCity of the channeC eCectrons in the 
moderate inversion region. ZJsing the proposed mo6ifity modeC, a new drain current modeC for the 
quicf^caCcuCation of the drain current in the moderate inversion region has aCso been deveCoped in 
this wor^. In addition to this, a new unified drain current modeC vaCid for aCC the three regions 
(i.e., the weah^ inversion, the moderate inversion, and the strong inversion) is presented in this 
wor^. 3 ? esuCts obtained from the proposed modeC have been verified with the experimentaC data 
reported in the Citerature for 0.26 jim and 0.09 jam channeC Cength CMOSCFIiTCs, which show a 
perfect match between the exp erimentaC data and the simuCated modeC characteristics. 

JLnother important issue for the Cow-power-anaCog-IC appCications is the proper modeCing 
and extraction of the device threshoCd voCtage parameter. In this worhj a new modeC for the 
threshoCd voCtage at the ‘ onset of the moderate inversion region aCong with its extraction procedure 
have been proposed. The deveCopment of this extraction method is based on two new integraC 
functions and I K proposed in this worh^ which are insensitive to the drain and source series 


resistances of the device and the experimental noise error. The results obtained from the proposed 
model are compared with the recently reported experimental data for 0.26 /am and 0.09 jam channel 
length devices and its reliability (i.e., short and narrow channel effects , and the parasitic series 
resistance effect on the proposed extraction method) is also tested and compared with the other 
existing methods by means of the JUtM-StPICT circuit simulator. The results reveal that the 
determined threshold voltage values always meet the ' onset of the moderate inversion region’ 
condition , and that, this extraction procedure is insensitive to the drain and source parasitic senes 
resistances and the noise introduced during measurements . 
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Chapter One 


Introduction 


Throughout the past four decades, both the productivity and the performance of 
microelectronics have advanced at exponential rates unmatched in technological history. 
Consequently, microelectronics has become the principal driver of the modem information 
revolution and the ubiquitous microchip has had a profound and pervasive impact on our daily 
lives - enabling such advances as microelectronic wrist watches, hearing aids, implantable 
cardiac pacemakers, pocket calculators, personal computers, wireless cellular telephones, 
optoelectronic-fiber networks, communication satellites, and the Internet. Propelled by the rapid 
advances made in the fabrication technology over the last four decades, the number of transistors 
per microchip has skyrocketed to about 100 million or more, while the cost of a chip has fallen 
constantly. The ever shrinking device dimensions, made possible by the improvements in the 
photolithography techniques, are essentially responsible for such a high packing density. 

Computer Aided Design (CAD) and simulation of electronic circuits are the main factors 
contributing to the success of VLSI. For accurate simulation of electronic circuits, the devices 
used in the circuits must be modeled accurately. The modeling of MOS transistors for CAD has 
been driven by the needs of digital circuit designers for many years. However, when these same 
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models are used for analog work, the outcome is very different: analog circuit behavior, as 
predicted through the use of such models, can have serious errors [1]. Thus, the rapid growth in 
low-power, low-voltage analog IC design imposes tighter requirements on the model accuracy. 

The continuous scaling down of MOS devices due to rapid advancements made in VLSI 
technologies require analog circuits to operate with a power supply voltage of less than 5 Volts, 
with a strong tendency to go further down to 3.3 Volts and below in the near future. At the same 
time, the development of high-performance portable equipment (e.g., radio receivers, disk drive 
filters, medical aids, etc.) requires the design of analog circuits having low power drain. The 
power consumption in analog ICs can be greatly reduced if the amplifiers are designed to work 
in the weak inversion region (as this region offers a low drain current) and, thus, the thermal 
cross talk between the devices gets reduced and consequently, the packing density gets boosted. 
However, the low drain current in this region also tends to degrade the speed of the circuit. 
Therefore, as the operating region of the device approaches the weak inversion, the packing 
density and the low power consumption will tradeoff with the speed of the circuit. Hence, m 
most analog circuits, the best tradeoff among area, power, and speed can be achieved when the 
devices work in the moderate inversion region, which lies in between the weak inversion region 
and the strong inversion region. With analog VLSI gaining popularity day by day, the moderate 
inversion region of operation of MOSFETs is becoming increasingly important for analog circuit 
designers. Therefore, an accurate modeling of the drain current in this region becomes an 
essential step in low power mixed-signal circuit design. 

Another important feature of the modem low power analog IC applications is that the 
devices are typically biased just above the threshold voltage (Vx) in order to reduce the dc power 
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dissipation. Thus, the threshold voltage is becoming a key parameter in MOSFET device design, 
characterization, modeling, and simulation, particularly for modem devices with very small 
geometry and relatively low power supply voltage. Hence, an accurate definition and extraction 
of this parameter is extremely important. 

1.1 Device Models 

Mathematically, a device model is an approximation of the device behavior, which can be 
measured or obtained from device characterization. It is a set of mathematical equations, which 
describe the physical behavior of the device under various input conditions. The accuracy of a 
model is determined from the deviation of the behavior predicted by the model from the 
measured behavior, e.g., from the current-voltage (I-V) characteristics. The reasons why 
MOSFET modeling for circuit simulation is still precarious are manifold. One reason is the 
improvement in the fabrication technology: device models based on long-channel MOSFETs 
with high power supply voltages are losing their validity currently, simply because modem high- 
performance MOSFETs feature extremely short channel lengths and are operated with a lower 
power supply voltage. Another reason is the interaction between the circuit complexity and the 
sensitivity to modeling errors: small errors in the measured device terminal characteristics and 
even smaller errors in their derivatives may turn into large errors during circuit simulation,, or 
even worse, these may create convergence problems in the simulations. In devices facricated 
from advanced deep-submicron CMOS processes, various short-channel effects occur, e.g., the 
Drain Induced Barrier Lowering (DIBL) effect, the dependence of the threshold voltage (Vr) on 
the channel length (L) and the channel width (W), etc. It is virtually impossible to account for all 
these effects accurately with a simple, physically based model. Consequently, the number of 
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model parameters increases, with most of these being merely fitting parameters without having 
much physical justifications. 

1.2 The Deficiency of the Existing Mobility Models 

The carrier mobility in the inversion layer of a MOSFET is an important parameter in the 
simulation of its drain current characteristics. The model for the effective carrier mobility in 
silicon inversion layers and in the bulk of the semiconductor is the one, which mostly affects the 
accuracy of the results of the terminal current calculations for a MOSFET. The effective mobility 
(p. e ff) of the inversion layer electrons is mainly dependent on the Coulomb scattering effect (i.e, 
the scatting due to the ionized impurities), the phonon scattering effect (i.e., the scattering due to 
the lattice vibration), and the surface roughness scattering effect (i.e., the scattering due to the 
surface states) [2,3]. It is found [2] that, in the moderate inversion region, jj, e ff is mainly 
determined by the Coulomb scattering effect, whereas in the strong inversion region, the phonon 
and the surface roughness scattering effects play the dominant role in order to govern the 
effective mobility (jj. e ff)- It is also found that the Coulomb scattering effect is a linear function of 
the inversion charge density [2,4], Therefore, p e ff increases with an increase in the gate to source 
voltage (Vgs) in the moderate inversion region, where the effects of the other two scattering 
mechanisms are not that marked. However, as the phonon and the surface roughness scattering 
effects are decreasing functions of Vgs [2], hence, p e ff decreases with increasing Vgs in the 
strong inversion region [5]. The BSIM mobility model [6], which is widely used for circuit 
simulation, considers only the phonon and the surface roughness scattering effects and ignores 
the effect of the Coulomb scattering. As a result, although the BSIM model gives quite accurate 
results in the strong inversion region, however, it seems to lose its validity in the moderate 
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inversion region. Hence, for low power analog IC design purpose, it is extremely important to 
define a mobility model, which is valid both in the moderate inversion region as well as in the 
strong inversion region. 

1.3 Parameter Extraction 

A device model is an abstract mathematical representation of the actual physical 
phenomena occurring in the device, which can be used to predict the behavior of the device 
under known stimulation with a reasonable accuracy. In order to apply a particular device model 
in circuit simulation, it is necessary to find a practical way of extracting the model parameters 
from the experimental data. The accuracy of the parameter extraction procedure depends on the 
accuracy of the model and on the accuracy and completeness of the experimental data. The 
parameter extraction procedure should be such that it gives fairly reasonable results without 
parameter adjustment, and can produce a nearly perfect fit to the device characteristics, and, 
hence, predict the circuit performance fairly accurately using popular circuit or device simulators 
(e.g., AIM-SPICE [7], MICROTEC [8], MEDICI [9,10], etc.). 

1.4 Challenges of the Threshold Voltage Extraction 

To date, several definitions [8,11] and methods [9,10,12-14] have been presented in the 
literature in order to define and extract the threshold voltage (Vt) parameter. The popular 
methods frequently used for V T extraction from the measured drain current (Id) versus gate-to- 
source voltage (Vgs) characteristics are 

1 ) the constant current (CC) method [12], 
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2) the transconductance change (TC) method [13], 

3) the ratio method [14], 

4) the difference function (D) method [10], and 

5) the transconductance (g ra ) linear-extrapolation (GMLE) method [9]. 

For designing low-power analog circuits, it is required that the definition of Vt and its extraction 
method should be 

1) reliably applicable to the moderate inversion region, i.e., the region where the carrier 
concentration in the inversion layer starts to become comparable to the depletion charge 
density, 

2) insensitive to the parasitic source/drain (S/D) series resistance and the surface roughness 
factor, 

3) insensitive to the measurement error due to noise, etc., and 

4) easily applicable to analog CAD tools. 

Unfortunately, however, none of the existing Vt extraction methods, as illustrated earlier, satisfy 
these requirements simultaneously. Most of these methods use the derivative method (i.e., taking 
the derivative of Id with respect to Vgs) and the linear extrapolation technique (i.e., extrapolating 
the linear portion of the g m versus Vqs characteristics) in order to extract Vj from the 
experimental data. However, both these techniques have some serious disadvantages. Firstly, if 
there is a measurement error due to noise, etc., then the error value will go higher after taking the 
derivative. Secondly, the extrapolation technique is highly sensitive to the S/D resistance and the 
surface roughness factor [13]. Since this technique depends upon the slope of the linear region of 
the Id versus Vgs (or the g m versus Vgs) characteristics, it is obvious that this slope changes with 
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a variation in the parasitic series resistances. As reported by Wong et al. [13], the slope of the 
linear region of the Id-Vgs characteristics also depends upon the scattering caused by the 
surface roughness, since it degrades the effective mobility of the channel electrons. Therefore, 
the extracted value of Vt will be dependent upon the value of the S/D resistance and the surface 
roughness factor. On the other hand, the D-Method [10], which uses the integral function, is free 
from the parasitic resistance effect and the measurement errors. Thirdly, most of these methods 
(except the GMLE method) use the drain current model, which does not take into account the 
Coulomb scattering effect in the moderate inversion region. As the Coulomb scattering effect is 
the main factor which determines the effective mobility of the carriers in the channel in the 
moderate inversion region, hence, these dram current models fail to explain the device behavior 
in this region accurately. Moreover, they only consider the above threshold region (i.e., the drift 
region) of the experimental data in order to extract the threshold voltage. As a result, the 
extracted value of the threshold voltage cannot properly represent ‘the onset of the moderate 
inversion’. Rather, the threshold voltage extracted from these methods represents the ‘onset of 
the strong inversion’. Therefore, these methods are not suitable for the threshold voltage 
extraction for devices to be used in low-power analog circuit designs. 

1.5 Goals 

The moderate inversion region of a MOSFET is a special region of operation [1], with 
properties distinct from those of the strong inversion and the weak inversion regions. Hence, this 
region should not be considered just as a “transition” region between the strong inversion and the 
weak inversion regions, which is mainly used to maintain the continuity of the two 
characteristics in the two regions. Modem low power analog IC design demands a more accurate 
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and robust model for the dram current m the moderate inversion region. The goal of the present 
work is to provide an accurate model of the effective inversion layer electron mobility, which is 
valid both in the moderate inversion region as well as in the strong inversion region. The existing 
BSIM mobility model is taken up in this work and modified accordingly in order to make it 
applicable for both the moderate and the strong inversion regimes. The results obtained from the 
proposed model is then compared with the recently reported experimental results for 0.26 pm 
[15] and 0.09 pm [16] channel length devices, and this comparison shows a perfect match 
between the two. 

The present work also addresses all the problems mentioned above with regard to the 
threshold voltage extraction for submicron channel length MOSFETs. In this work, an improved 
and unambiguous definition for the threshold voltage along with its accurate extraction 
procedure have been proposed, which takes into account the Coulomb scattering effect in the 
moderate inversion region, and, thus, is able to extract the value of the threshold voltage at the 
‘onset of the moderate inversion region’. The threshold voltage extraction method developed in 
this work uses both the drift and the diffusion components of the drain current in order to 
determine the threshold voltage. Moreover, it uses the integration function rather than the 
differentiation techniques adopted by the earlier authors, and, thus, is insensitive to the S/D 
series resistance effect, and suppresses the experimental errors due to noise, etc. The universality 
and reliability of this method is demonstrated by simulation studies using the AIM-SPICE 
simulator with varying device parameters (e.g., gate length (L), gate width (W), and S/D series 
resistances (R s and Ra)), as well as by comparing the results with the reported experimental data 
for 0.26 pm [15] and 0.09 pm [16] channel length devices. The results have shown a very good 
match among our simulation studies, the AIM-SPICE simulations, and the experimental data. 
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Chapter Two 


Theoretical Background 


The theoretical background required for an understanding of this work is briefly 
presented in this chapter. It highlights the basic concepts associated with MOSFETs, e.g., the 
threshold voltage model, the drain current model, the subthreshold region of operation, the 
scattering mechanisms, the mobility model, etc. 

2.1 The Electron Mobility 

Under low electric fields, the carrier drift velocity v is proportional to the electric field S. 

Thus, 

v = M‘- (2-1) 

The coefficient of proportionality po is called the low field carrier mobility. This equation can be 
obtained from the second law of motion for an electron moving in an electric field, given by [17] 

m n77 = c l^- m n-> ( 2 - 2 ) 

dt x 

where m* is the electron effective mass and q is the electronic charge. The first term in the right- 
hand side of Eq.(2.2) represents the electron acceleration by the electric field, and the second 
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term describes the effects of the collisions of the electrons with the lattice vibrations, impurities, 
and crystal imperfections on this motion. This term limits the electron drift velocity and the 
electron drift momentum. The term x appearing in Eq.(2.2) is called the momentum relaxation 
time, and, usually, it is of the order of 10 -12 to 10” 14 s. Under steady state, the electrons attain a 
constant velocity, given by 

v = qx-4- (2.3) 

m n 

Under low electric fields, x and m* are independent of the electric field, hence, the carrier 
velocity v can be given by v = po<= , where the low field electron mobility po can be expressed as 


Mo =q-V* 

m n 


(2.4) 


It should be noted that this linear dependence of the electron drift velocity on the electric field 
does not hold good under high electric fields, where the electrons may gain a considerable 
amount of energy from the electric field. This process can be described by the following equation 


dE 

— = aS • v ■ 
dt 


/ E-E 0 ^ 

V X E 


(2.5) 


Here E is the electron energy, Eo [= (3/2)kT] is the electron energy under the conditions of 
thermal equilibrium, and te is the effective energy relaxation time, and is usually of the order of 
10“ 11 to 10“ 13 s. The low field mobility po is determined by the electron collisions with phonons 
and impurities. These collision processes are called the scattering mechanisms, and are discussed 
in the next section. 
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2.2 The Scattering Mechanisms 


An electron moving in a periodic potential essentially behaves as a free electron, with the 
periodic potential only changing its mass from the rest mass to the effective mass. A perfectly 
periodic potential is, however, an idealistic assumption. This arrangement of the atoms in a 
perfect crystal is periodic, however, all the atoms move at all times about their equilibrium 
positions due to the thermal vibration at a given temperature. This vibration of the atoms changes 
the periodic nature of the potential. In addition, there are impurity atoms in semiconductors, 
either intentionally introduced (doping) or inadvertently got introduced during the various 
processing steps. The potential produced by these impurity atoms also alters the periodicity of 
the potential distribution. The potential distribution in a crystal can also be altered by crystal 
imperfections. Lastly, the description of the electron is built on the assumption that there is a 
single electron moving in the penodic potential. Other electrons when present are also expected 
to cause perturbation in the periodicity of the potential in their immediate vicinity. Due to these 
anomalies m the periodicity of the potential discussed above, an electron in a solid cannot 
execute long range motion without changing its state, which occurs when it arrives near a 
scattering center (e.g., impurity atoms, thermal vibrations of the lattice atoms, structural defects, 
etc.). The process causing the electron to change its state is called the collision process or the 
scattering process. 

In MOSFET inversion layers, the effective mobility (p e ff) of electrons mainly depends 
upon three scattering phenomena: p c , due to the Coulomb scattering; p p h, due to the phonon 
scattering; and p sr , due to the surface roughness scattering [2,3]. These three scattering effects 
are discussed as follows. 
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2.2.1 The Coulomb Scattering Effect 

There are different types of Coulomb scattering centers that can affect the carrier 
mobility in the inversion layer. The role of these Coulomb centers can be described as follows 

2.2.1. 1 The Ionized Impurity Scattering 

An ionized impurity center in the lattice of a semiconductor produces a long-range 
Coulomb field with a potential energy U, given by U = q 2 /(4nsr), where s is the permittivity of 
the medium and r is the radial distance between the electron and the Coulomb center. The motion 
of a carrier in such a field is in many respects similar to the analogous problem of Rutherford 
scattering of charged particles by ionized nuclei [18]. As shown by Takagi et al. [2], the effective 
mobility of electrons decreases with an increase in the substrate doping density. Hence, the 
ionized impurity scattering increases with an increasing doping density. 

2.2.1.2 The Remote Charge Scattering (RCS) 

The aggressive scaling down of MOSFETs is quickly approaching the deep submicron 
regime. For sub- 100 nm devices, an equivalent SiC >2 thickness of less than 2 nm will be required, 
according to the National Technological Roadmap for Semiconductors (NTRS) [19]. This is due 
to the fact that, as the dielectric constant of SiCh is low (« 3.9), one has to scale down the oxide 
thickness in order to boost the oxide capacitance value (since the conductance of devices in 
analog circuits is proportional to the oxide capacitance). However, the serious disadvantages of 
the ultrathin oxides are 
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1) Since SiC >2 is covalent in nature, the dopant atoms from the heavily doped poly gate can 
easily diffuse through the dielectric to the semiconductor layer, and, hence, can alter the 
doping concentration of the MOSFET channels. 

2) The gate leakage current (due to the direct tunneling of electron from the gate to the channel) 
can change the normal current voltage characteristics of the device [20]. 

Recently, alternative gate dielectrics (e.g., Ta20s, S 13 N 4 , AI 2 O 3 , etc., [19]) have demanded 
extensive attention, since they are expected to combine with (or even replace) SiC >2 in order to 
reduce the problems mentioned above. These alternative gate oxides have the following 
properties. 

• Since these materials have high values of permittivity (e), hence, a higher value of oxide 
capacitance per unit area (C' ox ) can be obtained with a relatively high value of oxide 

thickness (t ox ), and, obviously, the gate tunneling leakage current would reduce for a thicker 
dielectric. 

• Most of these alternative gate oxides are ionic in nature, and, hence, the dopant atoms can not 
easily diffuse through it. 

Therefore, as suggested by Momose et al. [20], the development of high-e gate insulators for 
MOSFETs is important as a future replacement of ultrathin SiC> 2 . Another disadvantage of 
ultrathin gate oxides is that the Coulomb scattering effect from remote charges in the poly gate 
(i.e., the ionized atoms which are introduced in the poly gate in order to make it conductive) or 
from poly gate/oxide interface roughness may lead to mobility degradation of channel electrons 
in MOSFETs having ultrathin gate oxides. The RCS effect on the effective mobility of the 
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inversion layer electrons may not be significant for higher oxide thicknesses (> 3 nm or so) or 
even for the alternative gate oxides, since the distance between the inversion layer electrons and 
the remote charges on the poly gate becomes too large for these cases to initiate the RCS effects. 

Besides the substrate impurities and the remote poly gate charges, the interface state 
charges and the charges trapped in the SiC >2 layer can also act as Coulomb centers, and take part 
in the Coulomb scattering event. It is shown that [2], when the substrate impurity concentration 
(Na) is greater than 5 x 10 16 cm -3 , the Coulomb scattering is dominated by the ionized impurity 
scattering and p c becomes inversely proportional to N A . However, when N A is less than 5 x 10 16 
cm -3 , other Coulomb centers determine the mobility degradation due to Coulomb scattering. 

2.2.2 The Phonon Scattering Effect 

The scattering of the inversion layer electrons by the vibration of the lattice atoms is 
known as the lattice scattering or the phonon scattering. In a semiconductor device, there exist 
two systems: one is the electron system and the other is the lattice system. The energy gained by 
the electron system from the applied electric field is transferred to the lattice system via the 
emission of phonons. There are two types of phonons: the acoustic phonon and the optical 
phonon. The dispersion characteristics of these phonons are shown in Fig.2.1. (The energy 
versus momentum characteristics of these phonons will be similar to their dispersion 
characteristics, since the former can be found by multiplying the ordinate and abscissa of the 
later characteristics with h and h respectively). From this figure, it can be noticed that the energy 
of the optical phonon is larger than that of the acoustic phonon, however, the variation of this 
energy (contrary to the acoustic phonons) is very little with the momentum change [21], Under 
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FIGURE 2.1 : Measured phonon spectra in silicon as a function of the wave 
vector along the [001] direction. The legends TO, LO, TA, and LA stand for 
transverse optical mode, longitudinal optical mode, transverse acoustic mode, 
and longitudinal acoustic mode respectively (after Sze [21]). 


15 




low electric fields, the emission of low energy acoustic phonons is responsible for lattice 
scattering and under large electric fields, the lattice scattering is due to the high energy optical 
phonons, since these optical phonons are not excited in thermal equilibrium due to the large 
amount of energy required for this process to happen. However, at very high electric fields, the 
carriers will gain enough energy between successive collisions in order to excite optical phonons. 
This process of emission of optical phonons sets an upper limit on the energy which a carrier can 
gain, since during optical phonon emission, a small change of energy requires a large change in 
the momentum. Because of this energy limitation, the drift velocity tends towards saturation, 
and, consequently, the effective mobility decreases [18,22]. However, this lattice scattering 
effect is more important for the heterojunction (e.g., GaAs-AlGaAs) MESFETs than Si- 
MOSFETs. 

2.2.3. The Surface Roughness Scattering Effect 

At the Si-Si02 interface, the periodicity of the lattice ceases to exist, and various foreign 
materials (which may come from contaminations during the fabrication process) may get 
attached to the surface. As a consequence, the number of intermediate levels within the forbidden 
band gap at the insulator-semiconductor interface will increase sharply, and the recombination 
rate can be greatly increased. In a practical MOS structure, there exist four types of surface states 
or charges, which can trap electrons from the inversion layer [21]. These are 

1) the surface states or the interface states, which are defined as the energy levels within the 
band gap at the Si-Si02 surface, and which can exchange charges with the semiconductor in a 
short time, 
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2) the fixed surface charges, which are located near or at the semiconductor surface and are 
immobile under applied electric fields, 

3) the mobile ions, such as sodium, potassium, etc., which are mobile within the insulator layer 
under the bias temperature stress conditions, and 

4) the ionized traps, which can be created, for example, by X-ray radiation. 

Since the number of carriers recombining at a surface will be dependent on the area 
rather than the volume, it is convenient to define a recombination rate per unit area in the 
following way 


r sr = sn s or r sr = sp s , (2.6) 

where r sr is expressed in carriers recombining per unit area per second, and s is called the surface 
recombination velocity (a measure of the “quality” of the surface) [18]. Equation (2.6) is 
equivalent to saying that the excess carrier concentrations nandp are moving towards the 
surface with a velocity s, i.e., there is a diffusion of excess carriers towards the surface where 
they disappear by recombination. Therefore, the higher rate of recombination at a surface 
reduces the average surface lifetime, which is anyway lower than the bulk lifetime. As a matter 
of fact, when the current flows near the surface, the scattering due to surface recombination can 
be extremely important. As the gate voltage increases, the electrons are attracted more and more 
towards the Si-SiC >2 interface, and the surface roughness scattering becomes more dominant. 


2.3 The Screening Effect 

The screening effect is an important phenomenon, which controls the mobility of the 
inversion layer electrons in the moderate inversion region. To understand this phenomenon, let 
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us consider a point charge q placed in a medium with permittivity 8 Then, the Coulomb potential 
cp 0 (r) at a distance r from the point charge can be expressed as 


<Po( r ) = 


q 

4irsr 


(2.7) 


Now, if the same point charge is placed in a sea of conduction electrons having a concentration 
of no, then the expression for the modified Coulomb potential cp(r) can be given by [23-25] 


cp(r) = -3— exp(-k s r), (2.8) 

4 Tier 

where k s , which is given by the relation k 3 = 4(3/7t) 1/3 ng /3 /a 0 (here ao is known as the Bhor 
radius), is known as the screening parameter and has the unit of cm -1 . The potential cp 0 (r) is 
known as the unscreened potential and cp(r) is known as the screened potential. The screened and 
the unscreened potentials, which are given by Eqs.(2.7) and (2.8) respectively, are schematically 
plotted as a function of the radial distance in Fig. 2.2. As can be seen from this figure, the 
screened Coulomb potential ( cp(r) ) behaves similar to the unscreened Coulomb potential 

(cpo(r) ) in the vicinity of the atom, but falls off drastically for large values of r due to the 
screening of the charge by the electron cloud. 


2.4 The Effective Mobility 

Figure 2.3 shows the dependence of the electron mobility in the inversion channel on the 
effective electric field (E e fi) at 300 K and 70 K for different values of the substrate doping 
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Potential 



Radial Distance 


FIGURE 2.2: Schematic diagram for the comparison between the screened 
and the unscreened Coulomb potential of a static charge, as a function of the 
radial distance. 
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densities (Na), as reported by Takagi et al. [2], where E e ff in the inversion channel can be written 
as [2] 

E eff =e;‘(xQ;+Qi), (2-9) 

where s S i is the permitivity of silicon, Q' is the inversion charge density per unit area, Q' d is the 

bulk (depletion) charge density per unit area, and x is a fitting parameter with typical values 
being 0.5 for (100) electrons and 0.33 for holes and electrons in (11 1) or (1 10) planes [2], These 
experimental results confirm that the electron mobility follows the universal relationship with 
Eeff (which is given by the relation p e ff E“° f 3 and p e ff « E~j? f at 300 K and 77 K respectively 

and are shown by dotted line in Fig.2.3) [2] for lower substrate doping densities, however, for 
higher substrate doping densities, it exhibits a significant deviation from the universal curve 
near the low field region or the moderate inversion region. The cause behind such deviation is 
believed to be due to the Coulomb scattering effect [2], which becomes more important when the 
doping concentration is high and the gate voltage or the normal field is low, since in the low field 
region, p c is much lower as compared to p S r or p P h. There is less effect of the Coulomb scattering 
on mobility when the inversion charge density is high because of the charge screening effect, 
which has been discussed in the earlier paragraphs. 

Figure 2.4 shows the dependence of the different mobility components on the effective 
field (E e ff). According to this figure, the effective mobility characteristics can be broadly divided 
into three contributing terms, i.e., the phonon scattering term (p P h), the Coulomb scattering term 
(p c ), and the surface scattering term (p sr ). According to Matthiesian’s rule [3], the effective 
nobility (p e ff) is usually assumed to be related in a reciprocal manner to the individual 
contributions by these three mechanisms as [2,3] 
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FIGURE 2.3: The effective electron mobility (p, e ff) at 300 K and 77 K in n-channel 
MOSFETs versus the effective normal field E e ff, for varying substrate acceptor 
concentrations (N A ), as given by Takagi et al. [2]. 
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Mobility 



FIGURE 2.4: Schematic diagram of the effective electric field (E e ff) 
dependence of the carrier mobility in the inversion layer by the three dominant 
scattering mechanisms, i.e., the phonon scattering (p P h), the surface roughness 
scattering (p sr ), and the Coulomb scattering (p c ), as given by Takagi et al. [2]. 
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1 


( 2 . 10 ) 


M-eff Me Mph Msr 

When E e ff is low, i.e., the channel is in the moderate inversion region, the value of (J. c is much 
lesser than either p. p h or p. sr . Hence, in this region, the effective mobility is determined mainly by 
the Coulomb scattering (i.e., p e f f « jj. c , since p c « p P h, Msr). A further study of this curve reveals 
that in this region, the Coulomb scattering term (p c ), and, subsequently, the effective mobility 
(Heff) increases with an increase in the effective electric field E e ff (i.e., with the gate voltage). In 
Fig.2.5, the mobility variation due to the Coulomb scattering effect has been plotted as a function 
of the inversion layer electron density (N') for different values of the substrate doping density 
(Na). From Fig.2.5, it can be seen that the mobility due to the Coulomb scattering term (|i c ) 
increases almost linearly with the inversion layer electron density (N') and decreases with 
substrate doping density (Na). This phenomenon is due to the screening effect of the electrons in 
the inversion layer [2,26], which has been discussed in Section 2.3. 

As stated previously, in a MOS structure, there are many charge centers located near the 
Si-Si0 2 interface, called the Coulomb charge centers, and they play important roles in carrier 
scattering. If the Coulomb centers near the interface are positively charged, then they attract 
electrons, and, in general, tie them up to a bound state. The properties of the bound states depend 
strongly on the amount of screening by the inversion layer electrons. At low temperatures, in the 
absence of screening, the bound states are deep enough to trap the first electrons that enter the 
inversion layer. As the gate voltage is increased, more and more electrons are added to the 
inversion layer, and some of these electrons will enter the quantum subband and will contribute 
to screening, thus weakening the scattering potential of the Coulomb centers. Alternatively, it 
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FIGURE 2.5: Inversion carrier density (N[) dependence of p c for n-channel 
MOSFETs with varying substrate impurity concentrations (Na), after Takagi et al. [2]. 




can be stated from Eq.(2.8) that, as the inversion layer electron concentration increases with 
increasing Vqs, the screening effect by the inversion layer electrons also increases. As a 
consequence, the scattering effect due to the Coulomb centers gets reduced and the effective 
mobility of the electron increases with the gate to source voltage (V os ) until the effect of the 
Coulomb scattering becomes comparable to that of the phonon scattering. With higher V GS (he., 
when the channel is in the strong inversion region and the screening is too high due to the high 
electron density), jj, c becomes very much larger than p. P h and p. sr . Hence, according to the 
Matthiesian’s rule, it can be inferred that for high values of Vgs, the effect of the Coulomb 
scattering on p. e ff becomes insignificant, and it is mainly determined by the phonon scattering 
(p P h) and the surface roughness scattering (p sr ) effects. 


2.5 The Threshold Voltage (V T ) 

When a voltage V G is applied to the gate of a MOSFET, part of it appears as a potential 
drop V ox across the oxide and the rest appears as a band bending V P S in silicon. It is 
straightforward to write that [7] 


Vo=v„+'r s =^-+T„ 


( 2 . 11 ) 


where Q' is the total charge per unit area induced in silicon, and C' ox is the oxide capacitance 
per unit area given by C^ x (= e ox / t ox ) for an oxide thickness to X , and s ox = £oSr(ox), where eo is the 

permittivity of free space and Sr(o X ) is the relative permittivity of SiC> 2 . There is a negative sign 
before Q' in Eq.(2.11) because the charge on the metal gate is always opposite in sign to the 
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charge in silicon, i.e., Q' is negative when V G is positive and vice versa. Now, if the work 

function difference is <|) ms (= <J> m - (|> s , where (j> m is the metal work function and <j) s is the 
semiconductor work function), and the equivalent oxide charge per unit area at the oxide-silicon 
interface is Q ' m , then the flat-band voltage (Vfb) can be given by 


V re =4> 


.Qi 

ms p t 


Introducing the Vfb term in Eq.(2.1 1), it can be rewritten as 


q; 


V g =V fb+ 'F s --^- 


( 2 . 12 ) 


(2.13) 


The termQ' has two components, one is the depletion charge density per unit area (Q' d ) and the 
other is the inversion charge density per unit area (Q'). These can be modeled by the following 
equations [7] 


Q; = a*. . and 

(2.14) 

Q' _ |2S st kTn f c q*F s /2kT 

(2.15) 

V N A 


Here Na is the acceptor doping density in the substrate, q is the electronic charge, n ; is the 
intrinsic carrier concentration for Si, k is the Boltzmann constant, T is the temperature (in K), 
and e S i is the permittivity of silicon. Equation (2.13) allows us to calculate the threshold gate 
voltage (V T ). It is assumed that at V G = V T , the surface potential OF,) is pinned at a value equal 
to <|>b, and there is no further change in with an increase of V G . At this surface potential, the 
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charge of free carriers induced at the insulator-semiconductor interface is still small compared to 
the charge in the depletion layer (i.e.,Q' d » Q' ). Therefore, using Eqs.(2.13) and (2.14), the 
expression for V T can be formulated as follows 


- V FB + <t> B + 


V2SsiqN A c|) I 

CL 


(2.16) 


Different theories have been proposed in the literature in order to model <|)b for an accurate 
extraction of Vt, and some of these are listed below. 


1. The most conventional definition states that the threshold voltage occurs when <|>b = 2<)>f, 
where 4 >f is the bulk potential, given by 


4> f = v .h ln 


v n . j 


(2.17) 


Here v t h is the thermal voltage, defined as v t h = (kT/q). At this point, the surface becomes as 
strongly n-type as the bulk is p-type 

2. According to Tsividis [1 1], the surface potential (¥ s ) actually increases with Vo beyond the 
value of 2<|)f, however, the rate of this increase is much slower than that in the subthreshold 
regime (i.e., for V F S < 2fc). According to this theory [1 1], the surface potential TL saturates at 
a value equal to 2<|)f + (4~5)vth. At this value of x f' s , the channel becomes strongly inverted. 
From this theory, (|)b can be modeled as follows 

*J>B « 24>f + C4~5)vth - (2.18) 
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3. Recently, a new theory has been proposed by Conde et al. [8]. They have considered both 
the inversion and the depletion charge densities in order to calculate the threshold voltage. 
According to their theory, <j>B can be formulated by the following equation 


<1) B = 2(j) F + v th In 


^240 

l v thC. 


(2.19) 


Here C, is a fitting parameter. With the help of a 2-D device simulator (MICROTEC [8]), it 
has been demonstrated that using C, = 10 yields the closest agreement of the results obtained 
from this theory with the majority of the threshold voltage data obtained from the various 
extraction methods [9,10,12-14] for a wide range of substrate doping densities. With C, = 10 
and using standard values of substrate doping densities give tj>B « 2<|>f + 0.98v t h. Therefore, it 
seems that this threshold voltage model gives an average value of the threshold voltage, 
which is somewhere in between the ‘onset of the moderate inversion’ and the ‘onset of the 
strong inversion’ regions, as per Tsividis [1 1]. 


2.6 The Drain Current Model 

To compute the drain current Id, it is convenient to adopt the ‘regional’ approach, where 
the various operating regions of MOSFETs are identified by the values of the gate to source 
voltage (Vos) or by the surface potential 0F S ), as discussed below. 

2.6.1. The Drain Current Model in the Subthreshold Region 

The subthreshold region, by definition, is the operating region of the MOSFET, where 
0 < vj/ $ < 2 <j) F . in this region, the drain current is entirely diffusive in nature (i.e., it is primarily 
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contributed by the diffusion of electrons from the source to the channel). This subthreshold 
component of the drain current (I SU b) can be modeled in the following way [11] 


T - Z 

sub T M'n 


yc;, 


L 2/L5<j> F + V SB - V( 


2 ^'>HF/v tha (V GS -V T )/ n v, h n - V DS / v 


rY th e 


(i-e - DS ,h ), 


( 2 . 20 ) 


bos N 

where y is the body effect coefficient, given by y = and r[ is an ideality factor, 


CL 


defined by r\ = 1 + 


2^1.5*, +V ffl 


. Here p n is the effective mobility of the electrons (in 


cm 2 /V-sec) in the subthreshold region, Vsb is the source to substrate voltage, and Vds is the drain 
to source voltage. Now for any V D s, I SU b can be simplified to 


!sub = l o ex P 


V nv th y 


( 2 . 21 ) 


where Iq, which is independent of Vgs, can be expressed as 


r W 

I 0 = ~-p 


yCL 


y 2 ^0 5<j) F / V th ^ 0 -VDS ^ V th ^ . 


L n 2Vl.5()) F + V SB -v th Ih 


(2.22) 


2.6.2 The Drain Current Model in the Strong Inversion Region 

The operating region of MOSFETs, where Vgs > V T is defined as the strong inversion 
region, where the surface potential is pinned at a value equal to 2<j) F + (4~5)vth. In this region, the 
drift process contributes primarily to the current conduction. This strong inversion component of 
the drain current (Ii nv ) can be modeled with the help of BSIM [6] in the following way 
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w [ y2 \ 

r ,nv =— C> eff (V gs -V t )V ds --^L for (V GS - V T ) > V DS , (2.23) 

L V 2 

and 

W 

I,nv =— C> eff (V GS -V T ) 2 (l + ^V DS ) for (V GS - V T ) < V DS • (2.24) 


Here p. e ff is the effective mobility of the carriers in the channel, given by [6] 


fW = 


Mo 

i + e( v Gs _v T ) ’ 


(2.25) 


where p.o is the low field mobility and 0 (in V -1 ) is the mobility degradation coefficient. The 
parameter A, (in V -1 ) is known as the channel length modulation parameter. Here the velocity 
saturation co-efficient (Ui), which was proposed in BSIM, is omitted for the sake of simplicity. 
Now, for low values of Vos and Vgs, one can neglect the V GS and the 0(Vqs - Vt) terms m 

Eqs.(2.23) and (2.25) respectively. Thus, under these conditions, Eq.(2.23) can be rewritten in 
the following way 


I,nv=K(V GS -V T )V DS , 


(2.26) 


where K (in A/V 2 ) is known as the device transconductance parameter and is given by 
W 

K=— C> 0 . 
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Chapter Three 


Development of the Mobility Model 


In this chapter, a new mobility model for short channel MOSFETs has been developed, 
which is based primarily on the BSIM mobility model. The Coulomb scattering effect at low 
electric field has been considered in this new model in order to apply it in the moderate inversion 
region. 


3.1 The Definition of the Threshold Voltage 

Today the world of IC devices has become broadly divided into two categories: one 
belonging to the analog devices and the other to the digital devices. The topologies, challenges, 
and the design techniques of these two types of devices are completely different. For most digital 
IC operations, the MOSFET devices act as electronic switches, which can either be in the ON 
state or in the OFF state. Therefore, for digital IC simulation and characterization, an accurate 
modeling of the moderate inversion region and the subthreshold region is not that important. 
These regions are only studied in order to measure the OFF state current of the digital IC 
, devices. Therefore, for digital IC design and simulation, the definition of the threshold voltage 
(V T ) should be such that it represents the ‘onset of the strong inversion region’ (i.e., V T s), such 
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that when Vgs becomes greater than Vts, the device is assumed to be strongly turned ON, 
otherwise it is OFF. 

On the other hand, as analog VLSI gains popularity, designers are trying to bias the 
devices in the moderate inversion region in order to reduce the power dissipation. An immediate 
consequence is that the threshold voltage is becoming the most important parameter for analog 
IC design, simulation, and characterization. In analog ICs, the MOSFET devices act as variable 
resistances (contrary to the electronic switches for digital ICs). Therefore, an accurate modeling 
of the moderate inversion region is extremely important for devices in analog IC applications. 
For a proper modeling of the moderate inversion region, the definition of the threshold voltage 
should be such that it represents the ‘onset of the moderate inversion region’ (i.e., Vtm). 

Figure 3.1 shows the schematic diagram of the surface potential V F S versus the gate to 
source voltage Vgs characteristic for a typical MOSFET device. As per Tsividis[ll], the 
characteristic can be divided into three major regions: 

1) The weak inversion region (for V L < Vgs < V T m; where V G s = V L when ML = In this 
region, W s increases from <j>F to 2<j>F + A<J> almost linearly with respect to Vgs- The value of 
A<() (which is discussed in detail in Chapter 6) is a few millivolts. 

2) The moderate inversion region (for Vtm < Vgs < Vts): In this region, % varies between 24 »f 
+ A(j> to 2<j> F + (4~5)v th in a nonlinear fashion with respect to V GS - 
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FIGURE 3.1 : Schematic diagram of the variation of the surface potential x i' s with 
respect to the gate-to-source voltage V GS for a typical MOSFET device, as presented 
by Tsividis [111. 


3) The strong inversion region (for V T s < Vos): In this region, 'f's becomes almost pegged at a 

value equal to 24 >f + (4~5)v t h and does not show any significant change with any further 
increase in Vqs. 

To date, several drain current models [6,27] have been proposed in the literature in order to 
explain the MOSFET drain current characteristics. In these models, the threshold voltage 
parameter (Vt) acts as the boundary between the strong inversion and the weak inversion 
regions. Then effort has been made in order to find proper models for the strong inversion and 
the weak inversion components of the drain current. Afterwards, these two components have 
been joined using different types of smoothing functions [27] or by some interpolation 
techniques [28] in order to obtain the overall drain current model. However, these techniques are 
quite inadequate to model the moderate inversion region accurately due to the following facts. 

1) These models use the value of the threshold voltage, which is neither Vtm nor Vts, but an 
average of the two, and they only rely on the smoothing functions for the modeling of the 
moderate inversion region. 

2) These models use the BSIM mobility model, which does not consider the Coulomb scattering 
effect, and, hence, is inapplicable in the moderate inversion region. 

Therefore, these drain current models are more suitable for modeling devices in digital IC 
applications rather than for the analog IC ones. For analog IC design, in order to accurately 
model the moderate inversion region, one has to not only take Vtm as the threshold voltage, but 
also use a modified expression for the mobility in order to account for the Coulomb scattering 
effect. In conclusion, it can be stated that for digital and analog IC designs, designers have to 
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switch between two different threshold voltage (V T s and V T m) and two different mobility 
models. 


3.2 Development of the Modified Drain Current Model 

The BSIM drain current model as discussed in Section 2.6 is based on the following two 
major approximations, which are not valid in the moderate inversion region. 

1) It is based on the simple charge control model for the inversion charge density Q[ [i.e., Q| = 
Cox(Vgs -V t )], which obviously is not valid in the moderate inversion region, since is 

not a constant in this region, and also that the inversion charge density is not so high that one 
can neglect the depletion charge density in comparison with it. 

2) The BSIM mobility model [as given by Eq.(2.26)] takes into account only the phonon 
scattering and the surface roughness scattering effects, however, it ignores the Coulomb 
scattering effect, which is the most important parameter for determining the effective 
mobility of electrons in the moderate inversion region. 

Due to these two major approximations, the BSIM drain current model gives accurate results 
only when the device is biased in the strong inversion region, and fails in the moderate inversion 
region. Therefore, one has to modify the BSIM drain current model in order to apply it in the 
moderate inversion region. In this work, an effort has been made in proposing a new mobility 
model, so that the BSIM drain current model can be applied in the moderate inversion region 
without changing its basic form [i.e., In = (W/L)C' 0X p e ff (Vgs - Vt)Vds], but only substituting 
the BSIM mobility model with the proposed mobility model and Vt with Vtm 
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3.3 Modeling of the Inversion Charge 


Generally, the inversion layer electron density N' is modeled as [7,9,1 1] 

n ;=— (v os -v t „). (3.i) 

However, as presented by Hauser [3], this linear relationship between N' and V G s does not quite 
hold good in the moderate inversion region and gives appreciable error, since in this region ¥<; is 
not a constant but varies slowly with V G s. In the weak inversion region, as 'Fg varies almost 
linearly with Vqs, the inversion charge density increases exponentially with increasing Vos- On 
the other hand, in the strong inversion region, as 4^ starts to saturate, the inversion charge 
density becomes a linear function of Vos- Therefore, it can be stated that in the moderate 
inversion region, the variation of N' with respect to V G s is neither exponential nor linear. In this 
work, a better formulation of N' in the moderate inversion region has been developed, which 
can be expressed as 


N ; =5— (V os -V t „) 1! , (3.2) 

q 

where £, (in V -0 5 ) is a fitting parameter. The origin of this equation is semi-empirical. Physically 
it is obvious that the value of oN' /5V GS in the moderate inversion region is greater than that in 

the strong inversion region, and smaller than that in the weak inversion region. One point must 
be noted here that, this model is valid only in the moderate inversion region, and in the strong 
inversion region N[ follows the linear model [as given by Eq.(3.1)]. 
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3.4 The Proposed Mobility Model 


From the discussion of the previous chapter, it can be stated that m the moderate 
inversion region, the effective mobility of the inversion layer electrons is mainly determined by 
the Coulomb scattering effect (p c ); and in the strong inversion region, it is primarily controlled 
by the phonon scattering effect (p P h) and the surface roughness scattering effect (p sr ). From the 
work of Takagi et al. [2] and Manzini [4], it can also be stated that |o. c varies linearly with the 
inversion layer electron density (N[) in the moderate inversion region. Hence, it can be written 
that 


Pc*N; or f4 c ^(V GS -V TM ) 15 (3.3) 

Now, if one uses the term |a. p h sr in order to represent the phonon scattering and the surface 
roughness scattering effects, then with the help of the BSIM mobility model, p p h sr can be 
expressed as 


(3 - 4) 

Now, the overall effective mobility (p e ff) model of the inversion layer electrons should be such 
that it satisfies the following conditions 


P'eff ' 


for 

v TM 

< V GS ^ ^TS > 


M-eff ' 

w M'phsr 

for 

V GS 

> V TS • 

(3.5) 
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In order to obtain the overall mobility model, a new function U(Vgs) is developed in this work, 
which is defined as 


U(V GS ) = z 


y+ 


m(v gs -v T m) is 

m + (v gs -v tm )>- 5 


(3.6) 


where M and y are two fitting parameters having units of V 1 5 and z is an another fitting 
parameter having the unit of' V " 1 5 . The function U(V GS ) has been plotted as a function of (V GS - 
Vtm) for four different values of M (0.05 V 15 , 0.087 V 15 , 0.1 V 15 , and 0.17 V 15 ) in Fig.3.2. 
Now, for proper modeling of p e ff, the value of M has to be taken such that the second term in the 
square bracket of Eq.(3.6) approaches a constant value for V G s greater than Vts- This is due to 
the fact that the Coulomb scattering effect only dominates in the moderate inversion region. 
Generally, the difference between Vts and Vtm (i.e., the region of the moderate inversion) is 
within 0.2 to 0.3 V. Therefore, typical values of M range from 0.089 (= 0.2 1 5 ) to 0.165 (= 0.3 1 5 ). 
The values of the other two fitting parameters (z and y) have to be chosen in such a manner that 
the value of the quantity z x (y + M) tends to unity for V G s greater than Vts- Thus, the effective 
mobility model proposed in this work can be given as 


M'eff, proposed 


HqU(Vgs) 

1 + 0(V GS - V TM ) 


(3.7) 


From Eqs.(3.6) and (3.7) it can be seen that, in the moderate inversion region (i.e., for Vtm < V G s 
c V T s), the variation of the [1 + 0(V GS - Vtm)] term is quite small as compared to the variation 
of the U(V G s) function, since typical values of 0 lie between 0.6 and 0.8 V -1 . Therefore, in the 
moderate inversion region, the effective mobility (p e ff) increases with increasing V G s. For V G s > 
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V TS , the function U(V G s) tends to attain a constant value of unity, and the proposed mobility 
model shows a similar behavior as the well-known BSIM mobility model. Hence, it can be stated 
that the proposed mobility model incorporates the Coulomb scattering effect in the moderate 
inversion region, which is absent in the BSIM mobility model, and, therefore, the proposed 
model seems to be more reliable, accurate, and robust than the BSIM mobility model for 
modeling devices to be used in low power analog IC designs. 


3.5 The Drain Current Model in the Moderate Inversion 
Region 

In this section, an approximate drain current model applicable to the moderate inversion 
region is proposed, which can be very helpful for quick calculations of the drain current in the 
said region. Here, two new functions f(VGs) and Ii(Vgs) have been introduced, which can be 

defined as 


f(V GS ) - C(V G s V tm ) , and 


(3.8) 


h(V GS ) = 


U(V GS ) 

i+e(v GS -v TM )’ 


(3.9) 


where £ is a new fitting parameter having a unit of V 05 . For a set of typical values of z, M, y, 
and 0 (which will be discussed in detail in the next chapter), both h(V G s) and f(V G s) functions 
have been plotted as a function of (Vgs - Vtm) in Fig.3.3 for three different values of C (°- 85 
V -05 1 15 V -05 , and 2 V -0 ' 5 ). From this figure, it can be inferred that by a proper choice of the 
value of C, the function h(V G s) can be made to approximate the function f(V GS ) in the moderate 
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h(V os ) 

f(Vos) 


<* 

y _ o x 7—0 5 / ' 

V TM = 0.41V 
V TS = 0.65 V 
z= 10.39 V -15 
y = 0.012 V 1 5 
M = 0.0865 V 1 5 
0 = 0.66 V -1 




Moderate 

Inversion 


Strong 

Inversion 





inversion region [i.e., for V TM < Vqs < Vts, h(Vos) — > f(Vos)]. Thus, in this region, the proposed 
mobility model [Eq.(3.7)] can be approximated by 

M'eff ~ ~ V TM ) ' (3.10) 

This approximate mobility model is valid only in the moderate inversion region and loses its 
validity in the strong inversion region. Putting this expression of peff in Eq.(2.23) and neglecting 
the V^s term for low values of Vqs, the approximate dram current equation in the moderate 
inversion region can be given by 

W 

I,„, = Y C >o«V os -V t „) is V ds ■ (3.11) 

Similarly, this drain current model is also valid only in the moderate inversion region (i.e., for 
Vtm < Vgs < Vjs), since the mobility model (as given by Eq.(3.10)) only considers the Coulomb 
scattering effect and ignores the phonon and the surface roughness scattering effects. Although 
this is a semi-empirical approximate model, it is very useful for the hand-calculation of drain 
current in the moderate inversion region. In the next chapter, the validity of this drain current 
model has been tested with the experimental data reported in the literature. 
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Chapter Four 


Model Verification and Discussion 


In this chapter, the proposed mobility model is verified with some recently reported 
measured data for short channel MOSFETs. 


4.1 The Modified Drain Current Model 

The total drain current (I D ) can be expressed as the summation of the subthreshold 

\ 

component ( T sub ) and the drift component ( I' str ). Hence, one can write 

i =r +r • (4.i) 

X D x sub “ A str \ / 

Here, the subthreshold component I' ub is defined by [6,27] 


where I SU b can be given by 


T T 

rt A sub sat 

■^sub “ 


•^sub ^sat 


x sub 


; I 0 exp 


l r i v th J 


(4.2) 


(4.3) 
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In Eq.(4.2), I S at is a new fitting parameter, which represents the saturated value of the 
subthreshold current in the strong inversion region, andfits value can be extracted from the 
experimental data of the device. It is defined in order to remove the higher order discontinuity in 
the drain current equation at V GS = V TM [6,27]. By replacing all the (V os - V T ) terms in 

Eqs.(2.23), (3.6), and (3.7) with a new function Vgtmx and neglecting the Vq S term for low 
values of Vos, one can obtain the expression for the drift component I[ tt as 


W 

^str — — C ox fieff, proposed V GTMX V DS , 


where the new function V GT mx is defined as [27] 


Vgtmx — s 


( Vqs-Vtm \ 

1 + e s 

l J 


(4.4) 


(4.5) 


where s is a fitting parameter having the unit of Volts. This function is defined m order to join 
the weak inversion and the strong inversion components of the drain current at Vjm, and also to 
ensure the continuity of the transconductance and its higher order derivatives. Thus, the overall 
drain current expression becomes 


IsaJo e 


nvih 


M-o 


isat + I 0 e 


nv* 


W 

J- +—' C' 

fc- L 1 + 0 V, 


GTMX 


^ I V-Vgtmx 
V M + V^ x/ 


V V 

v GTMX v DS 


(4.6) 


The second term in the right hand side of the above equation can be approximated by Eq.(3.11) 
in the moderate inversion region. 
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4.2 Parameter Extraction 


In order to apply the present model in circuit simulation, it is necessary to find a practical 
way of extracting the model parameters from the experimental data. The accuracy of the 
parameter extraction procedure depends on the accuracy of the model and on the accuracy and 
the completeness of the experimental data. Basically, there are two approaches for the extraction 
of the model parameters from the experimental data. In the first approach, special measurements 
are performed, from which the parameters can be extracted directly. This is termed as the 
analytical approach, and Tsividis [11] described such a procedure in detail. However, this 
approach becomes impractical when the number of parameters to be extracted increases. The 
second approach that is used for such cases is the numerical approach, which uses optimization 
techniques. In this approach, the model parameters are extracted from a given set of data using a 
nonlinear least-square optimizer, which curve fits the model to the experimental data using the 
controlling parameters as the variables. 

In this work here, the nonlinear optimization method is used to extract the model 
parameters from the experimental data. The well-known GNUPLOT software [29] is used to 
carry on the nonlinear optimization procedure. For the experimental data, the Id- Vos 
characteristics for a 0.26 pm channel length n-MOSFET have been taken from the work reported 
by Langevelde and Klaassen [15]. These characteristics are measured at a constant low value of 
Vos (50 mV), and V G s is incremented from 0 V to 1.5 V in steps of 10 mV. The device has the 
following parameters: gate length (L) of 0.26 pm, gate width (W) of 10 pm, and oxide thickness 
(to X ) of 4.5 nm. The threshold voltage at the onset of the moderate inversion (V T m) is extracted by 
the method proposed in Section 6.1 of this work and it is found to be equal to 0.41 V. Then, 
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Eq.(4.6) is fitted with the experimental data using the nonlinear optimization software. In 
Fig.4.1, the Id-Vqs characteristics obtained from the proposed model are compared with the 
experimental data and a perfect match is seen between the two. The extracted values of the 
model parameter are given along side in Fig.4.1. 


4.3 Comparison between the BSIM Mobility Model and the 
Proposed Mobility Model 

In order to make a comparison between the BSIM mobility model and the proposed 
mobility model, the following steps are performed: 

1) The drain current model proposed by Deshpande and Dutta [27] (which uses the traditional 
BSIM mobility model) is taken up and fitted with the experimental data for a short channel 
device. 

2) The drain current model proposed in this work [as given by Eq.(4.6)] is then taken up and 
fitted with the experimental data for the same device. 

3) Then, the values of the model parameters obtained from steps 1 and 2 are used in order to 
make a comparison between the two mobility models. 

The details of the procedure are explained as follows. 

For the experimental data, the I D -V G s characteristics for a 0.26 pm channel length n- 
MOSFET have been taken from the work reported by Langevelde and Klaassen [15]. The drain 
current model, as proposed by Deshpande and Dutta [27], can be expressed as 
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FIGURE 4.1 : Comparison of the I D versus V os characteristics obtained from the proposed model wi 
the experimental data for a 0.26 pm channel length MOSFET, as reported by Langevelde and Klaasss 
[15]. This comparison has been shown both in log scale (in left) and in linear scale (in right). 
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(4.7) 


_ Usub W 


I +1 

sat ■ L sub 


+ T ^oxM'eff.BSIM^GTX^DS’ 


where I SU b is defined by Eq.(2.21), and |_i e ff and Vgtx can be defined as 


Heff.BSIM 


Ho 

i+evctx 5 


(4.8) 


and 


^gtx — s ln| 


s~Vt 'N 


1 + e 


(4.9) 


where s (in V ') is the fitting parameter. Here, it should be noted that the term Vj is used instead 
of Vtm because the authors [27] did not use the term Vtm in their drain current model. Thus, the 
overall drain current model becomes 


I D = 


i s aJo e 




isat +I 0 e 


nv, h 


W 

■ + —C' 

l ox i+ev, 


Ho 


-V V 

v GTX v DS 


GTX 


(4.10) 


Now, this equation is fitted with the experimental data using the nonlinear optimization 
software. The extracted values of po, 6, s, and Vt are 432 cm 2 /V-sec, 0.66 V" 1 , 0.0463 V, and 
0.503 V respectively. Using these model parameters, the BSIM mobility model [as given by 
Eq.(4.8)] is plotted as a function of Vgs in Fig.4.2, which shows a continuous fall off of p eff with 
increasing Vas- 
in the same figure, using the values of the model parameters (which are listed along side 
in Fig.4.1), the proposed mobility model [which can be found by replacing V G s - Vtm terms with 
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BSIM mobility model 
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W= 10 pm 
V DS =50nP 
V T m = 0.41 ' 
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FIGURE 4.2 : Comparison between the traditional BSIM mobility model with the model proposed in this 
vork, based on the experimental data for a 0.26 pm channel length MOSFET, as reported by Langevelde 
uidKlaassen [151. 




Vgtmx i n Eqs.(3.6) and (3.7)] is also plotted as a function of Vgs- From this figure, it can be seen 
that the proposed mobility model (p e ff (1 proposed) increases with V GS for V GS < V TM + 0.3 V. 
However, for the same condition, p e ff,BsiM decreases with Vgs- This behavior proves that the 
proposed model can explain the Coulomb scattering phenomenon in the moderate inversion 
region, whereas the BSIM mobility model seems to be overestimating the mobility in this region. 
As the gate voltage increases beyond V TM + 0.3 V, p e ff, proposed decreases with V GS due to the 
effects of the phonon scattering and the surface roughness scattering, and |i e ff,BSiM also shows a 
similar behavior in this region. For high values of Vgs, both the proposed model and BSIM tend 
to converge with each other. However, for high Vos, the two curves do not exactly coincide with 
each other since the values of Vp and s are quite different for the two cases. Therefore, it can be 
concluded that in the low field region (i.e., for V T m < V G s < V TM + 0.3 V), the proposed model 
differs from the BSIM mobility model due to the inclusion of the Coulomb scattering effect. 
However, for high electric fields, the Coulomb scattering has less effect on the effective mobility 
of the inversion layer electrons and p. e ff is mostly determined by the phonon and the surface 
roughness scattering effects. As the proposed model considers all the three scattering effects (p c , 
Pph, and p. sr ) and BSIM only includes the phonon and the surface roughness scattering effects, 
hence, these two models differ in the moderate inversion region, however, show a similar 
behavior in the strong inversion region. 


4.4 Verification of the Proposed Drain Current Model 

In this section, the drain current model, which is given by Eq.(3.11), is verified with the 
recently published experimental data for short channel MOSFETs. In this work here, nonlinear 
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optimization method is used to extract the model parameters from the experimental data. For the 
experimental data, the same I D -V G s characteristics for the 0.26 pm channel length n-MOSFET 
[15] have been taken and the extracted value of Vtm is found to be equal to 0.41 V. Using Eqs. 
(3.11) and (4.1), one can express the overall drain current model for the weak and the moderate 
inversion regions as 

I = V 4-T' 
a D 1 sub ^ ^str 

Vqs-Vtm 

T T e nv * w 

= — 0 V GS -V TO + Y C °*^ V °™ X V DS ‘ ( 4 - H ) 

Isa t +I«,e ^ 

Here, Eq.(3.1 1) is used in order to model the 1^. component of Eq.(4.1). The region of the 

experimental data is taken up for values of Vos from 0 V to 0.6 V, since the drain current model 
as given by Eq.(4.1 1) is valid only in the weak and the moderate inversion regions. Afterwards, 
this equation is fitted with the experimental data of the said region using the nonlinear 
optimization software. In Fig.4.3, the Id-Vqs characteristics obtained from the proposed model is 
compared with the experimental data and it shows a perfect match between the two. The 
extracted values of the model parameters are given along side in Fig.4.3. 

The results predicted by the proposed model are also verified with the experimental data 

reported in the literature for sub-0.1 pm channel length devices. As another example, the 

experimental data obtain from the I D -V G s characteristics for a 0.09 pm channel length device has 

been taken from the work reported by Mii et al. [16]. The device has the following parameters: 

0 

gate length (L) of 0.09 pm, gate width (W) of 10 pm, and oxide thickness (tox) of 35 A. The 

*TT o Sil' j '<(■ „♦ -is 1 4 4 
*» 
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parameter Vjm is extracted using the method proposed in Section 6.1, and is found to be equal to 
0.25 V. The region of the experimental data is taken up for values of V GS from 0 V to 0.45 V for 
the comparison between the model and the experimental data. Then, the parameter extraction 
procedure is carried out to find the values for the model parameters in the same way as was done 
for the 0.26 pm channel length device. In Fig.4.4, the Id-Vqs characteristics obtained from the 
proposed model is compared with this experimental data [16] and it shows a perfect match 
between the two. The extracted values of the model parameters are given along side in Fig.4.4. 
From these results, it can be inferred that the proposed model can be reliably applied for the 
characterization of the moderate inversion region behavior of short channel MOSFETs. 
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V GS (V) 

FIGURE 4.4 : Comparison of the I D versus V GS characteristics obtained from the 
proposed drain current model [Eq.(4.11)] with the experimental data for a 0.09 pm 
channel length MOSFET, as reported by Mii et al. [16]. 




Chapter Five 


Development of the Threshold Voltage 
Extraction Method 


Parasitic source/drain series resistance is frequently significant in submicron MOSFETs. 
It represents a complication for the extraction of the device model parameters from the 
experimental I-V characteristics. Many methods have been proposed in the literature to perform 
the threshold voltage extraction even in the presence of series resistance [9,12-14], Most of these 
methods rely upon differentiating the measured I-V characteristics, which tends to worsen the 
experimental measurement errors that might be present in the data. Conde et al. [10] has 
proposed an alternative method, which is based on the use of the integral of the data in place of 
the differentiation technique. They introduced a new function called the ‘D function’ [10,30], 
which acts as a low pass filter, and, thus, tends to lessen the effects of the measurement errors 
and the S/D resistance effects on the parameter extraction procedure. Based on their method, a 
superior method has been developed in this work in order to extract the threshold voltage (Vtm) 
at the ‘onset of the moderate inversion region’ from the experimental Id-Vgs characteristics. 
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5.1 The ‘Difference Function’ D 

Let d and r be two generalized two-terminal devices connected in series. A terminal 
voltage V is applied to this series combination causing a current I to flow through the devices, 
which produces voltage drops Vd and V r across them respectively. Let the current be defined by 
the functions f(Vd), h(V r ), and g(V) of the voltages across the terminals of the two devices and 
the series combination respectively. Thus, 

I = f (V d ) = h(V r ) = g( V) . (5.1) 

The voltages can be expressed by their corresponding inverse functions of current as 

g- 1 (I) = f-'(I) + h- 1 (I). (5.2) 

Using geometric considerations, equivalent to integration by parts, the following equalities can 
be written as 


JldV a + 
0 

jv d dl = I 0 V d0 , 

0 

(5.3) 

V r0 

JldV t + 
0 

Io 

jV r dl = I 0 V r0 , and 

0 

(5.4) 

V 0 I ( 

JldV + 

jv dl = I 0 V 0 • 

(5.5) 


o o 


where V d0 , V^, and V 0 are the upper voltage integration limits corresponding to an upper current 
integration limit lo at a certain point in the I-V characteristics of the devices and the terminals. 
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Integrating the sum of the two voltages of the devices with respect to the current from 0 to I 0 
gives 


‘o, *0 

Jvdl=j 


V„dl + 



(5.6) 


Substitution of Eqs.(5.3), (5.4), and (5.5) into Eq.(5.6) produces a similar expression for the 
integrals of the current with respect to the three voltages 


V 0 V d0 

J I dV = Jl dV d + 
0 0 



(5.7) 


Finally, subtracting Eq.(5.7) from Eq.(5.6), one can define the ‘difference function’ D as 



v i0 \ i 0 

J*V d dl - JldV d + [v, 
Vo o J Vo 


dl- 


^rO 

JldV, 


(5.8) 


This function is valid not only for this case of two devices in senes, but also in general for any 
number of generalized devices connected in series [10]. 


5.2 The Properties of the 4 D Function’ 

The first and second terms on the right hand side of Eq.(5.8) are defined as the content 
and co-content of the device, respectively, and V and I are the device terminal voltage and 
current respectively. The integral function can also be considered as a measure of the 
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nonlinearity of the device, which would be zero for a device with only linear elements, since the 
content and the co-content are identical for such a device [10]. 

Any arbitrary two-terminal device with both linear and nonlinear elements has the 
following properties [6] 

a) The summation of the contents over all branches (i.e., all the parasitic resistances and other 
nonlinear components associated with the equivalent circuit of the device) is zero. 

b) The summation of the co-contents over all branches is zero. 

Therefore, the summation of the function D over all branches is zero. In other words, D 
eliminates the effect of the linear elements, such as series resistances, and can be used for 
extracting intrinsic parameters of semiconductor devices. 

5.3 The ‘G Function’ 

The D function has two integral terms. In order to reduce the complexity of the calculations, a 
new function G can be derived from the D function in the following manner 

lo v o 

D = jV dl - Jl dV (5.9) 

0 0 

By partial differentiation of the first integration term on the right - hand side of Eq.(5.9), one gets 
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Here g D is defined as g D - where V DS is the dram to source voltage and Id is the dram 
current of the device. 

Although these functions are especially defined for three terminal devices, it can be used for the 
parameter extraction for two terminal devices as well. These functions have been derived 
empirically in order to eliminate the effect of the parasitic series resistances on the MOSFET 
model parameter extraction procedure. These equations can be used as a tool for MOSFET 
parameter extraction. Although these functions are purely empirical in nature, however, it can be 
proved (both theoretically and with the help of AIM-SPICE simulator) that these functions can 
eliminate the effect of series resistances from the I-V characteristics of any nonlinear device. 

5.5 Case of a Constant Series Resistance 

Let a current I, flowing through a constant resistance R, produces a voltage drop V across 
it. If V x and I x are calculated for the resistance R, one gets 

1 

8d “r • 

Putting this value of gD in Eqs.(5.12) and (5.13), one gets 



(5.14) 
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(5 15) 



From the above result, it can be inferred that if the V x and I x functions are calculated from the 
Vgs versus Id characteristics of a MOSFET, then it will eliminate the series source/drain 
resistances, provided these are constants. This property of the V x and I x functions can also be 
proved from AIM-SPICE [7] simulation carried out in this work. The level-10 MOSFET model 
in AIM-SPICE is used for this simulation, and the following parameters are used: a threshold 
voltage of 0.7 V, a mask channel width of 8 pm, a mask channel length of 0.5 pm, a junction 
depth of 0.10 pm with a lateral extent of 0.075 pm, an oxide thickness of 10 nm, a substrate 
doping of 10 17 cm -3 , and a low field mobility of 331.5 cm 2 /V-s. The value of V DS is taken to be a 
constant at 50 mV and Vgs is incremented from 0 V to 1.5 V in steps of 10 mV. Four Id-V G s 
characteristics have been simulated for four different values of source/drain series resistances 
(R s , R<j), and these are shown in Fig.5. 1(a). The V x and I x functions are numerically calculated 
for each of the four simulated characteristics. Then the numerically calculated values of the I x 
function are plotted as a function of the V x function for all the four cases, as shown in Fig.5. 1(b). 
It can clearly be seen that all the four characteristics appear almost superimposed on each other, 
and, thus, it is conclusively proved that the effect of the series resistance is eliminated by these 
V x and I x functions. The method of the threshold voltage extraction using these V x and I x 
functions is presented in the next chapter. 
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X (A) I D (A) 



FIGURE 5.1(a) : The Id~V gs characteristics of an n-MOSFET for four different values of series 
source/drain resistances, simulated with the help of the AIM-SPICE simulator (Level -10) [7]. 



0.00 0.20 0.40 0.60 0.80 1.00 1.20 1.40 

V X (V) 

FIGURE 5.1(b) : The I x -V x characteristics obtained from Fig.5.1(a) for the same values of 
series resistances as in Fig.5.1(a). 
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Chapter Six 


The Threshold Voltage Extraction 


In this chapter, an alternative technique for the threshold voltage extraction is presented. 
It uses the V x (V G s,Id) and I x (V gs ,Id) functions, which are proposed in the previous chapter. 

6.1 Description of the Procedure 

For the experimental data, the Id-Vqs characteristics for a 0.26 pm channel length n- 
MOSFET have been taken from the work reported by Langevelde and Klaassen [15]. The 
characteristics are measured at a constant Vds (50 mV) and Vqs is incremented from 0 V to 1.5 

V in steps of 1 0 mV. The device has the following parameters: gate length (L) of 0.26 pm, gate 

o 

width (W) of 10 pm, and gate oxide thickness (t ox ) of 45 A. For each value of Vgs, the 
corresponding values of V x and I x have been calculated numerically using Eqs.(5.12) and (5.13). 
For example, if Id - 2 x 10~ 7 A at V G s = 0.5 V, then V x (V gs ,Id) and Ix(V gs ,Id) at Vqs = 0.5 V can 
be calculated in the following way 
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The integration terms have to be calculated numerically. In this work, Simpson’s rule is used for 
the numerical integration. In this way, for each value of Vgs> one can get unique values of V x and 
I x . The function I x (VgsJd) shows a negative peak at Vqs = 0.41 V, when it is plotted as a 
function of Vgs> as shown in Fig.6.1. This peak is denoted by Vtm m this work. From the same 
figure, it can be seen that a similar negative peak occurs when I x (V G s,Id) is plotted as a function 
of V x (Vgs,Id)- The value of V x (Vgs,Id), at which this peak occurs, is denoted by in this 
work. From Fig.6. 1 , it can be seen that = 0.368 V. 


6.2 The Mathematical Foundation* 

Replacing the term (W/L)C' 0X poC with a new parameter Ki (i.e., Ki = KQ, one can 
rewrite Eq.(4. 1 1 ) as 


___ 1 sat A 0 


Vgs-Vtm 

I T e 
in'-* 


Vgs-Vt, 


4KV 15 V 

T XVi V rjTKyfY v 


1 v GTMX v DS 


T +T e nVth 

i sat + i 0 C 


( 6 . 1 ) 


* For a complete derivation of the equations, please see the Appendix. 
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V GS ,V X (V) 


numerically from the experimental data, as reported by 
on ofV G s and V x . 





In order to obtain a simpler dram current model, the following assumptions have been made m 
this work. 

V GS~ V TM 

1) For Vos > V TM , I 0 e ’ lv,h » I*, and, hence, the first term of the right-hand side of Eq.(6.1) 
becomes approximately equal to I sat . 

2) For Vgs < Vtm, Vqtmx = 0, and otherwise it is equal to (Vgs - Vtm)- 
Therefore, one can simplify Eq.(6.1) into the following equations 


V GS~ V TM 



forV GS < V TM , and 

(6.2) 

=’I„,+K 1 (V 0S -V TM )‘ S V DS , 

for V GS < V TM < V TS • 

(6.3) 


With the help of Eqs.(6.2) and (6.3), one can derive the expressions for I X (V GS ,I D ) as follows 

1 x = 7TT cx P( bV cs )[! “ bV DS ] - , for V GS <V TM , (6.4) 

bV os bV DS 


where a = I 0 exp(- 


V 

v TM 

vy* 


) and b = 


1 


and for V G s > V T m, 


bV, 


DS 


V, 


-(V, 


GS 


'TM 


)-!„ + 


K,(V os -V m ) 


25 


' Kj (V GS — V TM ) l 5 V, 


DS 


2.5 


DS 


(6.5) 


In the same way, the analytical expression for V x for Vgs ^ Vtm can be found as 
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V X = V GS - T1V* 

1-exp 

( v W 

V GS 


_ 

l OvJ 


( 6 . 6 ) 


and for V G s »'nv t h, 


V x * V GS -riv th . (6.7) 

The exact analytical expression for V x , for Vqs > Vtm, is not important for this work (since this 
work primarily aims at determining the exact value of Vtm) s and, hence, it is not derived here. As 
stated in the previous section, a negative peak in I x appears at V GS = 0.41 V, when I x (V gs ,Id) is 
plotted as a function of V GS (Fig.6.1). Here, an assumption is made in this work that the value of 
V GS corresponding to this peak denotes the ‘onset of the moderate inversion region’ (i.e. V TM = 

0. 41. V). Taking this value to be that of Vtm> the parameter extraction procedure is carried out in 
order to find the values of the model parameters to be used in Eqs.(6.4) and (6.5). The parameter 
extraction procedure is performed in the following way. 

1. One region of the experimental Id-V G s data is taken for values of V GS between 0 V and 0.41 
V. Then, Eq.(6.2) is fitted (by nonlinear optimization software GNUPLOT [29]) with the 
experimental data for the said region. The following values of a, b, and r| have given the least 
error between the experimental data and the results predicted by our subthreshold drain 
current model: a = 1.3183xlO -10 A, b = 26.75 V -1 , and r| = 1.43. 

2. Again another region of the experimental I D -V GS data is taken where V GS is varying from 
0.41 V to 0.61 V (since Eq.(6.3) is valid only within this region). Then, Eq.(6.3) is fitted (by 
nonlinear optimization software [29]) with the experimental data for the said region. The 
following values of I sat and Ki have given the least error between the experimental data and 
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the results predicted by our dram current model: I sat = 6.615xl0' 6 A and Ki = 1.488xl<T 2 
AV“ 25 

Now, these values of the model parameters are plugged into Eqs.(6.4) and (6.5) and the values of 
I x are calculated and plotted as a function of Vqs, as shown in Fig.6.2. T his characteristic can be 
treated as the theoretical I x versus V GS characteristic, and in Fig.6.2 it is denoted by A. The 
Ix(VgsJd) versus Vcs characteristic, which is calculated numerically from the experimental data 
and appears in Fig.6.1, is also plotted in Fig.6.2 and is denoted by B. This characteristic can be 
treated as the experimental I x versus Vgs characteristic. In Fig.6.2, the numerically calculated 
Ix(V GS ,Id) characteristic is also plotted as a function of the numerically calculated V x (Vgs,Id) and 
is denoted by C. From Fig.6.2, it can be seen that curves A and B show an almost exact match 
between them everywhere except near the peak region, and both of them exhibit negative peaks 
at V GS =0.41 V. 


6.3 Analysis 

For typical values of t] (which is almost always less than 2) the product term r|v th (= b -1 ) 
is always less than Vos (50 mV). Consequently, the bV D s term in Eq.(6.4) is always greater than 
unity. Therefore, below the onset of the moderate inversion (i.e., V G s < Vtm), the theoretical I x 
[Eq.(6.4)] values are always negative and decreases exponentially with an increase in V G s- The 
experimental I x values (numerically calculated) also show the same behavior for values of Vgs 
below Vtm- However, from Eq.(6.5), one can see that beyond threshold (Vqs > Vtm), the 
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I X (A) 



FIGURE 6.2 : The I x function calculated numerically from the experimental data as reported by 
Langevelde and Klaassen [15] as a function of V GS and V x (curves B and C respectively). The 
results simulated from the theoretical model of I x [Eqs.(6.4) and (6.5)] proposed in this work is also 
shown as a function of V GS (curve A) for comparison. 
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theoretical values of I x increase with increasing values of V GS This is because the subthreshold 
component of the drain current tends to saturate after the threshold voltage (Vtm) and the drift 
component of the drain current dominates the current conduction process through the inversion 
channel. The experimental values of I x also increase with increasing Vgs beyond the peak. As a 
result, both the theoretical (A) and the experimental (B) characteristics exhibit peaks, which 
separate the two regions, i.e., the weak inversion region and the moderate inversion region. As it 
was previously assumed that Vtm = 0.41 V, and both the experimental and the theoretical values 
of I x exhibit peaks at V G s = 0.41 V, therefore, it can be inferred that the assumption is valid and 
the negative peak in the I x function (while plotted as a function of Vgs) always occurs at Vgs = 
Vtm (i*e., at the onset of the moderate inversion). 

Now from Fig.6.2, it can be seen that the theoretically calculated values of I x (curve A) 
match with the numerically calculated values of I x (curve B) everywhere except near the 
threshold (Vos = Vtm) region. Again from the characteristic C, it can be seen that a similar peak 
in l x (while plotted as a function of V x ) occurs at V x « 0.368 V. This value of V x is denoted by 
V 24 , in this work. It is interesting to note that the theoretical values of I x deviate from the 
experimental values of I x starting from V GS = V 2 $. In Fig.6.3, the results obtained from the 
subthreshold drain current model, which is given by Eq.(6.2), is plotted along with the 
experimental Id-Vgs characteristics, with the help of the extracted values of the model 
parameters. From this figure, it can be seen that the results obtained from our subthreshold drain 
current model match with the experimental data for V G s ^ V 2 *, however, they differ afterwards. 
This behavior can be explained in the following way. The sub threshold current model, given by 
Eq.( 6 . 2 ), can also be written in the following manner [1 1] 
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Id (A) 



V GS (V) 


FIGURE 6.3 : Comparison of the I D versus V GS characteristics obtained from the proposed 
subthreshold component of the drain current model [Eq.(6.2)] with the experimental data for 
a 0.26 pm channel length MOSFET, as reported by Langevelde and Klaassen [15]. These 
two characteristics exhibit a close match till V GS = V 2((l = 0.38 V, however, they differ 
afterwards. 
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(6.8) 


I D = 


( W' 




M C oxYV t 2 


2 xSow 


e rv s (v GB )- 24 > F ]/, 
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‘ V SB /V, h 


*~Vrm ! v i! 


), 


where V G b is the gate to substrate voltage, V SB is the source to substrate voltage, and V DB is the 
drain to substrate voltage. The surface potential (T / s ) increases almost linearly with Vgs until it 
reaches a value of 2 (|>f. After this, 'Fs increases very slowly with Vgs and saturates at a value of 
2(j) F + (4~5)v th (at V GS = V TS ). From Eq.( 6 . 8 ), it can be said that for W s < 2<j)F, Id increases 
exponentially with V GS - However, for > 2<j> F , increases very slowly with V GS , and, 
consequently, the subthreshold current loses its exponential nature and tends to saturate. Since 
the experimental and the theoretical values of I x match for values of V G s upto V 2(t , and differ 
afterwards, it can be hypothesized that at V GS = V 24 , the surface potential 0 F S ) reaches a value of 
24>f- Another important feature that can be observed from Fig.6.2 is that the difference between 
V T m and V 2 $ is nearly 42 mV. Now from Eq.(6.7), it can be stated that the difference between 
Vjm and V 2(t , is approximately equal to qv t h (37.3 mV). Hence, it can be written that 


Vtm w V 2 ^)+ T)V t h. 


(6.9) 


This relation gives a new model for the threshold voltage at the ‘onset of the moderate 
inversion’. Now, V 2(i) can be formulated by putting (j) B = 2<j> F in Eq.(2.16). Therefore, in this work, 
the proposed model of Vtm can be given as follows 


^TM ~ ^FB + + r , + ^l^th 


(6.10) 
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According to this definition, V TM is only 30 to 40 mV higher than V 2l j ) . From this model, it can be 
hypothesized that the ‘onset of the moderate inversion’ does not actually occur at *¥ s = 2<t> F , 
rather it occurs when C F S is a few tens of millivolts greater than 2<|> F . 

It has been stated earlier that the theoretical and the experimental values of I x do not match 
near around the threshold (V G s = V TM ) region. This happens because of the following facts. 

• When Vgs varies between V^ and Vtm> the subthreshold current seems to deviate from its 
exponential model and tends towards saturation. Therefore, Eqs.(6.2) and (6.4) lose their 
validity in this region. 

• Here it is assumed that the subthreshold current (I SU b) is saturated at a value equal to I sat for 
Vgs values beyond Vtm- In fact, when Vgs is varying between Vtm and Vts, 'Fs does not 
remain a constant but varies between 2<t> F and 2<j) F + (4~5)v t h. Consequently, I SU b also does 
not remain a constant at a value equal to I sa t but varies slowly with V G s in the moderate 
inversion region. However, this rate of variation of I SU b with respect to Vgs is too small, and 
when the drift component of the drain current becomes much larger than the sub threshold 
component, one can assume I SU b to take on a constant value. 

In Section 6.5, a new drain current model is proposed in order to match the results obtained from 
the theoretical I x function with the experimentally obtained values of I x . 
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6.4 Verification of the Proposed Method with the 
Experimental Results for Sub-0.1 pm Devices 

The results predicted by the proposed method are also verified with the experimental data 
reported m the literature for sub-0.1 pm channel length devices. For the experimental data, the 
Id-Vgs characteristics for a 0.09 pm channel length device have been taken from the work 
reported by Mii et al. [16]. These characteristics are measured at a constant V D s (50 mV) and 
V GS is incremented from 0 V to 0.5 V in steps of 10 mV. The V X (V GS ,I D ) and I X (V GS ,I D ) 
functions are calculated numerically from the experimental data in the same manner as discussed 
earlier. In Fig.6.4, the numerically calculated I x (V gs ,Id) function is plotted as a function of V GS 
(denoted by B) and also as a function of V x (V G s,Id) (denoted by C). The characteristic B shows a 
negative peak in I x at V G s = 0.25 V. Taking Vtm= 0.25 V, the parameter extraction procedure for 
obtaining the values of the model parameters to be used in Eqs.(6.4) and (6.5) has been 
performed in the same manner as stated previously. After the parameter extraction, the following 
values of the model parameters are found: a = 1.81 x 10~ 9 A, b = 26.05 V -1 , r\ = 1.47, Ki = 1.94 
x 10~ 3 AV" 2 ' 5 and I, at = 1.25 x 10“ 6 A. 

These values of the model parameters are plugged into Eqs.(6.4) and (6.5) and the 
computed values of I x (V G s,Id) are plotted as a function of V GS (denoted by A) in Fig.6.4. Both 
the theoretical and the experimental characteristics of I x (V G s,Id) [curves A and B] show peaks 
at V GS = 0.25 V. Also, curve C shows a peak at V x = 0.21 V. The difference between these two 
peaks is equal to 0.04 V, which is nearly equal to r|v th (0.038 V). From Fig.6.4, it can also be 
seen that the experimental I x (V G s,Id) versus V G s characteristic starts to deviate from the 
theoretical I x (V G s,In) versus V G s characteristic from V G s= 0.21 V. Therefore, the observations 
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V GS , V x (V) 


FIGURE 6.4 : The I x function calculated numerically from the experimental data as reported by 
Mii et al. [16] as a function of V GS and V x (curves B and C respectively). The results simulated from 
the theoretical model proposed in this work for I x [Eqs.(6.4) and (6.5)] are also shown as a 
function of V GS (curve A) for comparison. 


75 



made for this device are quite similar to that of the 0.26 pm device, given earlier. Hence, it can 
be stated that the proposed method is also valid for sub-0.1 pm devices. 


6.5 Further Analysis* 

Near the peak region, the discrepancies between the values of the I x function, calculated 
from the theoretical model and from the experimental data, can be reduced using the drain 
current model as given by Eq.(6.1). The expression for I x , calculated from this drain current 
model, can be written as 


— sat 

~ bV us 


In 


JU+a exp(bV GS ) 


^sat a 



■V. 


Y\ 


TM 




dV GS 


I sat x a exp(bV GS ) R 
Isat + a exp(bV GS ) 


sin 


1 + exp 


V - V 

v GS v TM 


Y\ 

v TM 

S )j 


1 5 


V, 


DS 


( 6 . 11 ) 


However, it is not possible to get an analytical expression for the integration term appearing in 
Eq.(6.11). Hence, in this work, this integration has been performed numerically. For the 
experimental data, the Id-V G s characteristics for a 0.26 pm channel length device have been 
taken from the work reported by Langevelde and Klaassen [15]. The value of the parameter V TM 
is found numerically in the same manner as discussed in Section 6.1, which yielded Vtm = 0-41 
V, and the parameter extraction procedure is carried out to find the values of the model 
parameters to be used in Eq.(6.1). The parameter extraction procedure is performed in the same 
way as stated in Section 4.4. Now, these values of the model parameters (which are listed along 

* For a complete derivation of the equations, please see the Appendix 
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side in Fig.4.3) are plugged into Eq.(6.11) and the values of I x are calculated and plotted as a 
function of V G s, as shown in Fig.6.5. The I X (V GS ,I D ) versus V GS characteristic, which is 
calculated numerically from the experimental data is also plotted m the same figure. From this 
figure, it can be seen that the values obtained from the model of I x match almost exactly with the 
values of I x calculated numerically from the experimental data. Both the theoretical model and 
the numerical analysis show negative peaks at V GS = 0.41 V. 

The results predicted by the proposed method are also verified with the experimental data 
reported in the literature for sub-0.1 pm channel length devices. For the experimental data, the 
Id-V G s characteristics for a 0.09 pm channel length device have been taken for V G s values 
between 0 V and 0.5 V, from the work reported by Mii et al. [16]. The value of the parameter 
Vtm is extracted numerically as stated in Section 6.1, which yielded Vtm = 0.25 V, and the 
parameter extraction procedure is carried out to find the values of the model parameters in the 
same way as was done for the 0.26 pm channel length device. The values of the model 
parameters are listed along side in Fig.4.4. Now, using these values of the model parameters, the 
I x function is calculated and plotted as a function of V G s, as shown in Fig.6.6. The I x (V G s,Id) 
versus V G s characteristic, calculated numerically from the experimental data is also plotted in the 
same figure. From this figure, it can be seen that the values obtained from the model of I x match 
almost exactly with the values of I x calculated numerically from the experimental data. In this 
case, both the characteristics show negative peaks at V G s = 0.25 V. From these results, it can be 
strongly concluded that the negative peaks in the I x -V G s characteristics always occur at V G s = 
Vtm (i-e. at the onset of the moderate inversion region). 
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Experimental I x 
Theoretical L 


L = 0.26 pm 
W =10 pm 


FIGURE 6.5 : The I x function calculated numerically from the experimental data as 
reported by Langevelde and Klassen [15] as a function of V GS • The results simulated 
from the theoretical model of I x [Eq.(6.11)] proposed in this work are also shown as a 
function of V GS for comparison. 








6.6 Extraction of V TS 

When V G s becomes equal to V TS , the device is at the ‘onset of the strong inversion’. The 
V TS extraction method developed in this work is termed as the ‘g m , max ’ method. In this method, 
the values of g m are calculated from the experimental Id-Vqs data by differentiating Id with 
respect to Vos numerically Now, the value of Vqs at which the value of g m attains a maxima , 
represents the onset of strong inversion (i.e., at this point, Vgs = Vts). From the discussion in 
Section 2.4, it can be stated that for Vqs < Vxs, the effective mobility (p e ff) increases with Vgs 
due to the Coulomb scattering effect. However, for V G s > Vts, the effective mobility (p eff ) 
decreases with Vqs due to the phonon and the surface roughness scattering effects. As a result, 
gm increases with Vgs below Vts, peaks at Vqs = Vts, and then decreases with Vgs 

This method of V T s extraction may contradict with other methods like the g m , max/3 [7] or 
other linear extrapolation methods. However, it seems that the other methods extract the value of 
the threshold voltage, which is neither Vtm nor Vts, rather an intermediate value between the 
two. This is due to the fact that the other methods do not consider the Coulomb scattering effect, 
which is dominant in the moderate inversion region. 

It is true that the differentiation method adopted in this work for the extraction of Vts 
may depend on the series resistance effect and the errors due to noise introduced during 
measurements. However, for low power analog IC operation, Vtm (onset of the moderate 
inversion) is much more important than Vts (onset of the strong inversion). Hence, the accuracy 
in the extraction of Vxs is not that much important for analog IC design. Therefore, in this work, 
much more importance is given on the extraction of an accurate value for Vtm rather than Vts. 
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Moreover, an accurate modeling of Vis seems to be a difficult task, since the assumption that 
Qa » Q! does not quite hold in this region. 

In the next chapter, the proposed method is verified with the simulated data for devices 
having various gate lengths and widths. The effect of the parasitic series resistance on the 
proposed method has also been observed. The performance of the proposed method has been 
compared with other existing threshold voltage extraction methods. 
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Chapter Seven 


Reliability Test of the Proposed Method 


The threshold voltage of short channel MOSFETs is dependent on the channel length (L) 
and the channel width (W) of the device. The S/D series resistances also have an important role 
to play in determining the drain current characteristics for short channel MOSFETs. In this 
chapter, all these short channel effects are discussed. Also, the proposed method of Vtm 
extraction is tested for devices having different channel lengths and widths. The results are 
compared with other threshold voltage extraction methods. The effect of the source/drain 
resistance on various existing threshold voltage extraction methods is also presented in this 
chapter. 

7.1 Effect of Channel Length on the Threshold Voltage (Vtm) 

MOSFETs with short channel lengths experience a shift in the threshold voltage with a 
change in the channel length (e.g., when L decreases, Vtm also decreases). An explanation for 
this shift in the threshold voltage for short channel devices can be given in terms of the charge 
sharing (or the charge conservation) concept [7]. This model [7] considers the sharing of the 
depletion charge in the regions where the gate depletion zone overlaps with the depletion zones 
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of the source and the dram contacts. This shared charge is balanced by the counter charges 
distributed between the gate electrode and the source and the drain contacts. As a result of this 
charge sharing, the effective depletion charge, which is induced by the gate voltage, is reduced. 
Hence, the model of Q d , as given by Eq.(2.14), no longer holds good. Including the charge 
sharing effect, Eq.(2. 1 6) can be rewritten in the following way [1 1] 


V TM ~ V FB + + W^B + V SB 


1 - 


a,8 

__ 




B + ^SB 


(7.1) 


where 8 is given by 


8 = 



, and V TM 


is the effective threshold voltage. Here ai is a 


dimensionless fitting parameter. With a decrease in the channel length (L), the threshold voltage 
also decreases. This phenomenon should be reflected in the extracted value of the threshold 
voltage from the experimental Id-Vqs characteristics. In order to examine the effect of the 
channel length on the threshold voltage extraction procedure, AIM-SPICE simulation has been 
performed on a representative device. The level- 17 MOSFET model of AIM-SPICE is used in 
the simulation, and the following parameters are used: a mask channel width of 5 pm, an oxide 
thickness of 10 nm, a substrate doping of 10 16 cm -3 , a channel doping concentration of 10 17 cm' 3 , 
and a threshold voltage of 0.7 V at zero substrate bias. For the other parameters (e.g., the low 
field mobility, the mobility degradation coefficient, etc.), the default values of the level-17 model 
are taken. The value of V DS is taken to be a constant at 50 mV and V G s is incremented from 0 V 
to 1.5 V in steps of 10 mV. Five I D -V G s characteristics have been simulated for five different 
values of channel length (e.g., 0.5 pm, 0.4 pm, 0.3 pm, 0.2 pm, and 0.1 pm). From these five 
simulated characteristics, threshold voltages have been extracted using the GMLE method [9], 
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the TC method [13], and the proposed method. The first two methods can be discussed as 
follows. 

1. The GMLE Method [9]: The simulated Id versus Vos characteristics are numerically 
differentiated with respect to Vgs in order to obtain the transconductance (g m ) versus Vgs 
characteristics. The linear region of the g m versus Vgs characteristics is extrapolated and its 
intersection with the V G s axis gives the threshold voltage. 

2. The TC Method [ 13] : The simulated Id versus Vgs characteristics are numerically 
differentiated with respect to Vos in order to obtain the g m versus V G s characteristics. Again, 
these characteristics are differentiated with respect to V G s in order to get the second 
derivative of g in with respect to Vgs (i-e., 9 2 g m / <9 Vq S ) versus Vgs characteristics. The value 

of Vgs, at which the value of d 2 g m /SV 2 s is maximum, is taken to be the threshold voltage 
of the device. 

In Fig.7. 1 , the extracted values of the threshold voltage from each of these methods are plotted as 
a function of the channel length. It can be seen that the results of the method proposed m this 
work are consistent with the theory, and also the extracted values of the threshold voltage using 
the proposed method are comparable to others, and lie intermediate between the values extracted 

by the other two methods. 

7.2 Effect of Channel Width on the Threshold Voltage (V TM ) 

MOSFETs with narrow channel widths also experience a shift in the threshold voltage 
with a change in the channel width (e.g., when W decreases, V T m increases). For aMOSFET, the 
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depletion region undei the gate is not limited to just the area directly below the thm oxide, rather 
some of the field lines emanating from the gate charges terminate on the ionized acceptor atoms 
on either side ol the gate along the width of the channel. These field lines constitute what is 
called a fringing field. If W is large, the part of the depletion region on the two sides of the gate 
along the width of the channel is a small percentage of the total depletion volume, and, thus, it 
can be neglected. However, with values of W of only a few micrometers or less, the side parts of 
the depletion region become a significant fraction of the total depletion volume. Thus, the total 
effective depletion charge under the gate is increased, and, therefore, the threshold voltage 
increases with decreasing channel width. Including this narrow channel effect, Eq.(2.16) can be 
rewritten in the following way [11] 


V'JM ~ V,. B +(j) B +y7<t>B + ^SB 


1 + 


a 2 8n 

2W 


V4>b+ V sb 


(7.2) 


where a 2 is a dimensionless fitting parameter. In order to examine the effect of the channel width 
on the threshold voltage extraction procedure, AIM-SPICE simulation has been performed on a 
representative device. The level-17 MOSFET model in SPICE is used for this simulation, and the 
same values for the various parameters are used as stated in the preceding section. Here, the 
channel length is kept constant at 0.5 pm, and four I d -Vds characteristics have been simulated for 
four different values of gate width (e.g., 10 pm, 8 pm, 6 pm, and 5 pm). From these four 
simulated characteristics, the threshold voltages have been extracted using the GMLE method 
[9], the TC method [13], and the method proposed in this work. In Fig.7.2, the extracted values 
of the threshold voltage from each method are plotted as a function of the channel width. It can 
be seen that the results of the method proposed in this work are consistent with the theory, and 
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FIGURE 7.1 : Simulated values of the threshold voltage extracted by the GMLE method 
[9], the TC' method [13], and the method proposed in this work as a function of the channel 
length (based on the results from AIM-SPICE simulation). 



FIGURE 7.2 : Simulated values of the threshold voltage extracted by the GMLE method 
[01 the TC method [13], and the method proposed in this work as a function of the channel 

width (based on the results from AIM-SPICE simulation). 
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also the extracted values of the threshold voltage using the proposed method are comparable to 
others, and lie intermediate between the values extracted by the other two methods. 

7.3 Effects of Source/Drain Series Resistance on the 
Extraction of the Threshold Voltage (V TM ) 

The channel ol a MOSFET is in series with two ‘parasitic’ resistances, one associated 
with the source and the other associated with the drain. A schematic diagram of the current-flow 
pattern in the source/drain region of a MOSFET is shown in Fig.7.3. The symbols (e.g., l c , S, p c , 
Xj, x c , and R S( j) that appear in Fig.7.3 will be defined in the subsequent sections. The total 
source/drain resistance can be divided into several parts: R ac is the accumulation layer resistance 
in the gate-source (or the gate-drain) overlap region where the current mainly stays near the 
surface; R sp is associated with current spreading from the surface layer uniformly across the 
depth of the source-drain junction; R S h is the sheet resistance of the source/drain region where the 
current Hows uniformly, and Rc 0 is the contact resistance (including the spreading resistance in 
silicon under the contact) in the region where the current flows into a metal line. 

Once the current flows into the metal line, there is very little additional resistance 
introduced, since the resistivity of aluminum (which is commonly used as the interconnect metal) 
is very low (p A i. « 3xl0~ 6 Q-cm [31]). In VLSI interconnects, the thickness of the aluminum 
interconnects is typically between 0.5— 1.0 pm. Therefore, the sheet resistivity is of the order of 
0.05 0/1 1. This is negligible in comparison with the channel sheet resistivity p C h, which is 
typically between 2000-7000 O/D, except when a long, thin wire is used as the interconnect. 
Figure 7.3 only shows the series resistance on one side of the device. The total source/drain 
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resistance per device is, of course, twice of that shown m F lg . 7.3. The accumulation layer 
resistance R ac depends on the gate voltage. Since it is not easily separable from the active 
channel resistance, R ac is considered to be a part of the resistance of the channel having an 

effective length (L dT ) of the device. The other components of the source/drain resistance can be 
modeled as discussed below [31], 


7.3.1 The Spreading Resistance (R sp ) 

The spreading resistance R sp is given by [31] 


R s . = ^-ln 0.75— 
SP TtW X. 


(7.3) 


were W is the device width, and Xj and x c are the junction depth and the inversion (or 
accumulation) layer thickness respectively. As shown in Fig.7.4, the current spreading takes 
place in a uniformly doped medium with resistivity pj. 


7.3.2 The Sheet Resistance (R sh ) 

The sheet resistance R s t, is given by [31] 

(7.4) 

where S is the spacing between the gate edge and the contact edge, and Rsd is the sheet resistivity 
of the source-drain diffusion, typically of the order of 50-500 Q/Q [31]. Since the equivalent 
value of Rsd is much smaller than that of p ch , this term is usually negligible if S is kept to a 
minimum, limited by the overlay tolerance between the contact and the gate lithographic levels. 
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FIGURE 7.3 : A schematic cross section showing the pattern of current flow from a 
MOSFET channel through the source or drain region to the aluminum contact [31]. The 
diagram identifies various contributions to the series resistance. 


i 

x, 


Resistivity pj 

FIGURE 7.4 : A schematic diagram showing the resistance component associated with 
the injection region where the ament spreads from a thin surface layer into a uniformly 
doped source or drain region [31]. 
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7.3.3 The Contact Resistance (R co ) 

The contact resistance R eo is given by [31] 


R 


CO 



(7.5) 


where l c is the width o( the contact window (Fig.7.3), and p c is the interfacial contact resistivity 
(in Q-cm 2 ) of the ohmic contact between the metal and silicon. 

Now, the effective gate voltage (Vos.eff) and the effective dram voltage (VDs,eff) can be defined as 
X.s.cti V <iS -I n R s and V DSeff = V DS - I D (R D +R S ), (7.6) 


where Rd and Rs are the drain and source parasitic series resistances respectively and are given 
by Rs (Rn) R«„ * R M > + R sh . Thus, these parasitic resistances degrade both the Vos.eff and V D s,eff- 
As a result, there is a degradation of the drain current for short channel devices with low intrinsic 
resistance. This effect is more severe for ultra-short channel MOSFETs. As the device sizes are 
scaled down, the effective source/drain diffusion cross-sectional area also goes down, and, 
consequently, the source/drain series resistance increases. Moreover, in ultra short channel 
MOSFETs, LDD (lowly doped drain) structure is used in order to eliminate the hot electron 
hazards. This LDD structure also tends to increase the parasitic series resistance. The effect of 
series resistances is shown with the help of AIM-SPICE simulation in Fig.5.1(a). How these 
characteristics are obtained is discussed in Chapter 5. From these characteristics, the threshold 
voltage is extracted using the GMLE method [9], the TC method [13], and the method proposed 
in this work. Then, for each method, the error in the extracted values for Vt is calculated using 

the following equation 
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V T ~v 

% error = • — j^ Rs *° T R d^ s =Q 


T lR d ,R s =o 


This crroi value is plotted in Fig.7.5 as a function of source/drain resistance for each of the 
extraction methods. Since the proposed method can filter out the parasitic resistance effect, 
hence, it shows the least error among the three methods, and, therefore, it can be stated that this 
method is more reliable for ultra-short-channel length devices where the series resistances are 

too high. 


7.4 Mathematical Complexity of the Proposed Method 

It is more time consuming to perform a numerical integration than a numerical 
differentiation on a digital processor. Therefore, the proposed method for the threshold voltage 
extraction may take longer time than the extraction methods based on differentiation techniques, 
when it is executed on a computer. However, the computational time of the proposed method can 
be greatly reduced if one properly chooses the integration limits of the proposed V x and I x 
functions. If V !ow and V hlgh are the lower and the upper limits of integration, then in this work it 
is found that V U)W = V TM - 0.1 V and V hlg h = V TM + 0.05 V is a sufficient condition for getting a 
negative peak at the onset of the moderate inversion region. The function I x is calculated from 
the experimental data for the 0.26 pm channel length device [15] for two cases: m the first case, 
Viow and Vh,gh are taken to be equal to 0 V and 1.5 V respectively, and is plotted in Fig.7.6 
(denoted by curve A). In the second case, Vi 0 w and Vhigh are taken to be equal to 0.31 V (— Vtm 
0.1 V) and 0.46 V (= Vtm + 0.05 V) respectively, and is also plotted in Fig.7.6 (denoted by curve 
B). Both the curves show negative peaks at V G s = 0.41 V. 
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Percent Error (%) 



FIGURE 7.5 : Comparison of the precent error in the extraction of the threshold 
voltages using the GMLE method [9], the TC method [13], and the method 
proposed in this work as a function of S/D resistance, based on the results of AIM- 
SPICE simulation. 
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The same procedure has been repeated with the experimental data for a 0.09 pm channel 
length device [16]. In this case, curve A in Fig.7.7 is found by taking Vi ow = 0 V and Vh, g h = 0.5 
V, and curve B in the same figure is found by taking Vi ow = 0.15 V and V hlg h = 0.3 V. It is to be 
noted here also that both the characteristics show negative peaks at Vgs = 0.25 V. Hence, it can 
be concluded that by a proper choice of the integration limits, the mathematical complexity and 
the computational time of the proposed method can be greatly reduced. 
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FIGURE 7.6 : The I x function calculated numerically from the experimental data, as 
reported by Langevelde and Klaassen [15], as a function of V GS for two different 
integration limits. 



V GS (V) 

FIGURE 7.7 : The I x function calculated numerically from the experimental data, as 
reported by Mii et al. [16], as a function of V GS for two different integration limits. 
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Chapter Eight 


Summary and Conclusion 


With the growing interest in mixed analog/digital circuit design in recent years, the use of 
MOS transistors has not only been restricted to digital circuit design alone, but also it has been 
extended to analog circuit design as well. With analog VLSI gaining popularity day by day, the 
moderate inversion region (i.e., the transition region between the subthreshold region and the 
strong inversion region) of operation of MOSFETs is becoming increasingly important for circuit 
designers. This region becomes particularly important with the continuous downscaling of MOS 
transistors towards submicron dimensions and the subsequent decrease in the power supply 
voltage. A point will soon be reached when the power dissipation alone will be the most 
important factor governing the performance, reliability, and the overall system cost. Therefore, 
this growth of low power analog circuit design has given rise to the need for developing accurate 
MOSFET models, which are capable of explaining the device behavior in the moderate inversion 
region. 

To date, several drain current models have been proposed in the literature in order to 
explain the MOSFET dc characteristics. Although these models are quite adequate to explain the 
device characteristics in the strong and the weak inversion regions, however, most of them have 
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explain the device behavior accurately in the moderate inversion region because of the 
following facts. 


These dram current models use the value of the threshold voltage, which does not represent 
the onset of the moderate inversion region. 

These models use different types of smoothing functions or interpolation techniques in order 
to join the weak inversion and the strong inversion components of the drain current near the 
threshold region, and these models depend on these so-called smoothing functions for 
explaining the moderate inversion region behavior. As a result, these models fail to represent 
the moderate inversion region characteristics accurately. 


• These models use the well known BSIM mobility model, which is quite accurate in the 
strong inversion region, however, it seems to lose its validity in the moderate inversion 
region. 

In consequence of these drawbacks, these drain current models are applicable only for the 
digital IC design purpose but not for the analog one. Also, an accurate model for the effective 
mobility of the inversion layer electrons in the moderate inversion region is extremely important 
for the design of low power analog ICs. 

Another issue, which is tremendously important for low power analog IC design, is the 
proper extraction and modeling of the threshold voltage parameter. The existing threshold 
voltage extraction methods suffer from the following shortcomings. 
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They can not extract the threshold voltage at the ‘onset of the moderate inversion region’, 

i.e., when the inversion layer electron density become equivalent to the substrate doping 
density. 

• These extraction procedures are highly influenced by the parasitic source-drain resistance 
effect and the experimental noise error. 

In this work, an effort has been made in order to propose a new improved mobility model and 
a threshold voltage extraction procedure in order to portray the moderate inversion region more 
accurately. A brief summary of this work is given below. 


8.1 Summary of the work 

As presented in literature [2,5,9], the Coulomb scattering effect is the most important in 
influencing the effective mobility (p e ff) of the channel electrons in the moderate inversion region. 
Due to this effect, p e ff increases with increasing values of Vos, in the moderate inversion region. 
This phenomenon is opposite to the results found from the widely used BSIM mobility model 
[6]. This is due to the fact that the BSIM mobility model does not consider the Coulomb 
scattering effect in the moderate inversion region, rather it considers only the phonon scattering 
and the surface roughness scattering effects in the strong inversion region. In this work, an 
attempt has been made in order to propose a new mobility model so that the device behavior in 
the moderate inversion region can be explained accurately. A new function U(Vqs) has been 
proposed in this work in such a manner that when the existing BSIM mobility model is 
multiplied with it then one gets the modified mobility model, which can explain the Coulomb 
scattering phenomena in the moderate inversion region and has no effect in the strong inversion 
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region. An effort has also been made in order to propose a new dram current model for the 
moderate inversion region. This model can be immensely helpful to the low power analog circuit 
designers for the hand calculation of the dram current in the moderate inversion region. 

Results obtained from the proposed model have been verified with the experimental data for 
a 0.26 pm [15] and a 0.09 pm [16] channel length MOSFET, which shows a perfect match 
between the two. Hence, it can be concluded that the mobility model developed in this work here 
is more reliable and accurate in describing the drain current characteristics in the moderate 
inversion region for short channel MOSFETs. 

In this work, a new threshold voltage extraction method has been developed. The 
difference (D) function, proposed by Conde et al. [10], has been modified in order to apply it for 
three terminal devices. Based on their work, two new functions V x (Id,Vgs) and I x (Id,Vgs) have 
been proposed in this work. From the experimental Id-Vgs characteristics of any MOSFET, the 
V x (Id,Vgs) and I x (Id,Vgs) functions are calculated numerically. A negative peak appears in the 
I x (Id,Vgs) function when it is plotted as a function of Vgs- It has been proved in this work that 
the gate to source voltage (Vgs) at which this peak occurs, represents the ‘onset of the moderate 
inversion’ (i.e. Vtm)- A similar negative peak in the I x (Id,Vgs) function appears when it is plotted 
as a function of V x (I D ,V G s). In this work, it is found that the separation between these two peaks 
in the I x function is approximately equal to r|v th Volts. From theoretical calculation, it is 
observed that this second peak appears when the surface potential 0F S ) becomes twice of the 
bulk potential (<t> F )- From this result, a new model for the threshold voltage (Vtm) have been 
proposed in this work, which can be expressed as 
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v ~ v j_ oa i -^4s SI N A (j)p 

V TM ~ V FB + 2 <P F + + T1V* 


CL 


The proposed threshold voltage extraction method is different from the other methods in a 
sense that it extracts the threshold voltage from the experimental Id-Vgs data at the ‘onset of the 
moderate inversion’. Therefore, the proposed method and model for the threshold voltage 
extraction is much more accurate and reliable for device characterization for low power analog 
IC design. 

Also, an attempt has been made to propose an alternative method for the extraction of Vts 
(i.e., the threshold voltage at the onset of the strong inversion). However, we have not put much 
effort for the accuracy of this Vts extraction procedure, since this parameter is not that important 
for low power analog IC design. 

In summary, it can be stated that the proposed method for the threshold voltage extraction 
has the following advantages over the other existing extraction methods [9,12-14]. 

• It extracts the threshold voltage at the onset of the moderate inversion region. 

• The extraction procedure uses the integral function, which acts as a low pass filter, and 
suppresses the error due to noise that may be introduced during the measurements. 

. The extracted value of the threshold voltage is independent of the source/drain parasitic 
series resistance. 

The results predicted by the proposed method is compared with the recently reported 
experimental results for 0.26 pm [15] and 0.09 pm <16] channel length devices. Also, its 
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performance has been studied and compared with the other methods with the help of AIM- 
SPICE simulator. The short channel and the narrow channel effects and the parasitic series 
resistance effect on the threshold voltage extraction method are also observed using the AIM- 
SPICE simulator. Although the proposed method has a computational complexity comparable to 
that of the linear extrapolation method, it has the advantage of being less sensitive to the 
measurement error and noise and being independent of the parasitic resistance effect. 

Scope for Improvement 

• The subthreshold slope model as proposed by Thakare and Dutta [32] is based on the 
assumption that the subthreshold slope parameter is independent of the variation of the gate 
to source voltage (Vgs). However, in reality, due to the Coulomb scattering effect, the 
subthreshold slope becomes a function of Vgs- Hence, there is a scope for improvement of 
that subthreshold slope model [32] using the results of the present work. 

• The proposed mobility model can be further improved by considering the transconductance 
overshoot effect for ultra short channel MOSFETs, with the help of the correction factor F v 
as proposed by Kolhatkar and Dutta [33]. 

• In this work, the numerical integration for the proposed threshold voltage extraction method 
h as been performed using the Simpson’s rule of the numerical integration. The computation 
time of the proposed threshold voltage extraction method can be highly reduced if any 
numerical integration method can be developed which is faster than the commonly used (e.g., 
Simpson’s rule, trapezoidal rule, etc.) integration techniques. 
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Appendix 


Here, an effort has been made in order to find the analytical expressions for the I x and V x 
functions, where use has been made of the drain current model [as given by Eqs.(6.2) and (6.3)]. 


Al. Calculation of I x in the Subthreshold Region 


From Eq.(6.2), it can be written that 


( V -V 
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where a = I n exp — , and b = 

l nv th J T i v t h 



A2. Calculation of I x in the Moderate Inversion Region 

Again, from Eq.(6.3), it can be written that 


I D = I sat + K, (V GS - V TM ) L5 V DS for V GS > V TM ■ (A.2) 


In order to calculate I x in the moderate inversion region, it is divided into two components: I x i and I X 2 , 
where I x i is the subthreshold component and I X 2 is the drift current component. The component I x i can 
be calculated in the following manner, noting that for Vos < Vtm, Isub = a e bv ° s , and for Vqs > Vtm, Isub 
~ Isat- 


V GS 

T xl = Jsd^V 0 S — I D 


TM 

J 


- V 

0 V DS 


V GS t 

-*-^V os + }^dV os -I„ b 


v 

v™ V DS 


[ a exp(bV G s ) dv + filLdV -I 
0 J V DS GS . J V n5 


V TM 

L 


bV 


-[exp(bV TM )-l] + -f-[V GS -V TM ]-I 


DS 


V, 


sat 


DS 


bV, 


1 [I,-a] + i-[V 0 S -V TO ]-I 


DS 


V, 


sat 


(A.3) 


DS 


The component I X 2 can be calculated in the following manner, noting that for Vgs < Vtm, hnv ~ 0, and 
for Vgs > Vtm, hnv = Ki(Vgs - Vtm) L 5 Vds- 
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Hence, the total current I x can be given by 
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A3. Calculation of V x at the Subthreshold Region 

The function V x can be calculated in the following manner. 
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A4. Calculation of I x from the Accurate Drain Current Model 

In this section, I x is calculated from the drain current model as given by Eq.(6. 1). Thus, 
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